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Abstract 

Following Jacobi's geometrization of Lagrange's least action principle, 
trajectories of classical mechanics can be characterized as geodesies on 
the configuration space M with respect to a suitable metric which is the 
conformal modification of the kinematic metric by the factor {U + h), 
where U and h is the potential function and total energy, respectively. 
In the special case of 3-body motions with zero angular momentum, the 
global geometry of such trajectories can be reduced to that of their moduli 
curves, which record the change of size and shape, in the moduli space 
of oriented m-triangles, whose kinematic metric is, in fact, a Riemannian 
cone over the shape space M* ~ S'^(l/2). 

In this paper, we show that the moduli curve of such a motion is 
uniquely determined by its shape curve (which only records the change of 
shape) in the case of /i 7^ 0, while in the special case of ft = it is uniquely 
determined up to scaling. Thus, the study of the global geometry of such 
motions can be further reduced to that of the shape curves, which are 
time-parametrized curves on the 2-sphere characterized by a third order 
ODE (cf. Theorem 3.9). Moreover, these curves have two remarkable 
properties, namely the uniqueness of parametrization and the monotonic- 
ity, as stated in Theorem 4.6 and Theorem 5.8, that constitute a solid 
foundation for a systematic study of their global geometry and naturally 
lead to the formulation of some pertinent problems such as those briefiy 
discussed in §7. 
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1 Introduction 



1.1 Local and global characterization of 3-body trajecto- 
ries 

The classical 3-body problem of celestial mechanics studies the local and global 
geometry of the trajectories of a 3-body system under the influence of the gravi- 
tational forces, or equivalently a conservative system with potential energy —U, 
where 

is the Newtonian potential function. Thus, when the particles have position 
vectors SLi = (xi,yi,Zi) with respect to an inertial frame, the trajectories are 
locally characterized by Newtons equation 

dU mimj , , miirik , , r- • 71 r-, ^ n-, 

miai = — = — 3-^(aj -ai)H 3 — (a^ - aj), {i, j,k} = {1,2,3} (2) 

where = |aj— ai| are the mutual distances, and (mi,m2,m3), rrij > 0, 
is the given mass distribution, assumed to be normalized so that X^mj = 1. 
Since the above equation is of order two, a trajectory is completely determined 
by the initial position and velocity of the particles - often referred to as the 
deterministic doctrine of classical mechanics. 
We use the following notation 

7 = 5^mi|a/, T=^^mi\ai\\ 12 = ^ m,(a, x a,) (3) 

for the (polar) moment of inertia, kinetic energy and angular momentum, re- 
spectively. These are the basic kinematic quantities, and their interactions with 
the potential function U play a major role in the dynamics of the 3-body prob- 
lem. For example, it is easy to deduce the classical conservation laws from the 
system (2), namely the conservation of total energy 

h = T-U, (4) 

linear momentum ^TOjaj and angular momentum ft. As usual, the invariance 
of linear momentum allows us to choose the inertial reference frame with the 
origin at the center of mass. Moreover, by differentiation of / twice with respect 
to time t and using (2) we get 

BU 

7 = 4T + 2 ^ a, • — = 4T - 2C/ = 2(C/ + 2h) (5) 
c/aj 

where we have used the fact that U is homogeneous of degree —1 as a function 
of the vectors a^. This is the Lagrange- J acobi equation. 

On the other hand, trajectories can also be determined as solutions of a 
suitable boundary value problem, and the simplest and most basic one is, for 
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example 



"For a given pair of points P, Q, what are those trajectories 7(f), to <t<ti, 

with 7(to) = P and 7(ii) = Q ?" 

Then solutions are found by applying an appropriate least action principle, 

which characterizes solutions as extremals of an action integral J{'j) among 
those virtual motions ^{t) with the given pair of end points, together with some 
additional constraints. 

Here we shall focus attention on the two least action principles due to La- 
grange and Hamilton, which are quite different but dual to each other : 

Lagrange : ^1(7) = / Tdt, fixed energy h (6) 

Hamilton : 72(7) = / {T + U)dt, fixed time interval [to,ti] (7) 

The motions t — > 7(t) are regarded as parametrized curves in the Euclidean 
configuration space 

M = {(ai,a2,a3);5]]rniai = 0} (8) 

Our aim, however, is to reduce the study of 3-body trajectories to a study of 
associated curves in a lower dimensional space, namely the interior configuration 
space, i.e., the moduli space 

M ~ R^/SO{3) w c (9) 

With the appropriate assumptions, one expects that the least action principles, 
as well as Newton's differential equation, can be pushed down to the level of 
M. In fact, one of our major results is that the study of 3-body motions with 
vanishing angular momentum further reduces to the analysis of specific curves 
on the shape space, which is the sphere 5^ C R'^. 

The final step of our program is, of course, the lifting procedure from the 
moduli curve ^{t) back to the actual trajectory j{t). But this is a purely ge- 
ometric construction which is well imdcrstood (cf. e.g. [5] or the following 
subsection) and it will not be a topic here. Briefly, the curve in M determines 
the trajectory in M uniquely up to a global congruence. 

1.2 Riemannian geometrization and reduction 

Classical mechanics up to present time is largely based upon developments re- 
lated to Hamilton's least action principle, involving Hamiltonian systems and 
canonical transformations. Geometrically speaking, the underlying structure is 
the symplectic geometry of the phase space. However, in this paper we shall 
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rather focus on the Riemannian geometric approach, based upon Jacobi's refor- 
mulation of Lagrange's least action principle. In his famous lectures [6], Jacobi 
introduced the concept of a kinematic metric 

ds^ = 2Tdt^ 

on the configuration space M, in terms of the kinetic energy T of the mechanical 
system. For example, in the case of an n-body system with total mass ^ mi = 1 

ds^ = 2Tdt^ = mi{dx^ + dyf + dzf) (10) 

i 

Now, for a system with kinematic metric ds^, potential function U and a given 
constant total energy h, set 

M^UM) = {p(^M;h + U{p)>0} (11) 
dsl ^ {h + U)ds^ 

where ds^ is the associated dynamical metric. Then by writing 

dsh = Vh + Uds = Vrds = V2Tdt 

Jacobi transformed Lagrange's action integral (6) into an arc-length integral 




and hence, in one stroke, the least action principle becomes the following simple 
geometric statement : 

" Trajectories with total energy h are exactly those geodesic curves in the 
space with the dynamical metric ds\ " 

Nowadays, such metric spaces are called Riemannian manifolds, and the 
global geometric study of geodesic curves is often referred to as the Morse theory 
of geodesies. In particular, we note that the dynamical metric ds^ is a conformal 
modification of the underlying kinematic metric ds^. 

In this geometric setting, the notion of "congruence class" is defined by 
the action of the rotation group SO{3), fixing the center of gravity (= origin). 
It acts isometrically on the configuration space {M,ds^) with the kinematic 
metric, and also on the modified metric space {M (^jj ^-^ , ds\) for any 50(3)- 
invariant potential function U . The corresponding <S'0(3)-orbit spaces inherit 
the structure of a (stratified) Riemannian manifold with the induced orbital 
distance metric, which we denote by 

(M, ds^), {M^u,h),dsl ), dsl = {h + U)ds\ (13) 

and similar to (11), for h negative the geodesies must stay inside the Hill's 
region, namely the proper subset 

Mh = M^u,h) = {p&M;h + U{p)>Q} (14) 
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By definition, the projection map 

TT : M ^ M = M/SO{3) 

is a (stratified) Riemannian submersion, where the horizontal tangent vectors 
at p € M are those perpendicular to the S'0(3)-orbit. They are mapped, by 
the tangent map dw, isometrically to the tangent space of M at p = 7r(p). 
Via the map tt there is a 1-1 correspondence between curves 7 in M and their 
horizontal lifting 7 in M (rcsp. M(jjj^-^), up to congruence. In fact, for a (virtual) 
motion ^{t), the property of being horizontal is equivalent to the vanishing of 
the angular momentum vector fl. 

On the other hand, the above metric ds^ on M also has a kinematic inter- 
pretation in analogy with (10), namely 

ds^ = 2fdt'^ , f = T - T'^ (15) 

where is the purely rotational energy and hence the difference T, representing 
that of the change of size and shape, is naturally the kinetic energy at the level 
of M. Therefore, we also refer to ds^ (rcsp. ds\) as the kinematic (resp. 
dynamical) metric on M. Classical mechanics, indeed, tells us how the term 
T*^ can be calculated from ft via the socalled inertia operator of the system; in 
particular, it follows that T"^ = if and only if $7 = 0. 

Now, assume = and let Uhe a, nonnegative and S'0(3)-invariant function 
on M. Then it is not difficult to see that both action principles (6), (7) can be 
pushed down to M. In the first case, using Jacobi's reformulation (12), we arrive 
immediately at the following geometric statement similar to the one above : 

" Curves in M representing 3-body trajectories with total energy h 
(and n = 0) are exactly those geodesic curves in the moduli 
space Mfjj^h) with the induced dynamical metric ds^." 

In the case of (7), the Lagrange function L = T + [/ is also defined at the 
level of M . Indeed, when $7 = we can view M as the configuration space 
for a simple conservative classical mechanical system, namely with potential 
energy — [/, kinetic energy T, and conserved total energy h = T — U. It is easy 
to calculate the associated Euler-Lagrangc equations with respect to suitable 
coordinates in M, as demonstrated in Section 3.2. Finally, the reduced Newton's 
differential equation on M can be calculated by the procedure described in [5], 
but we leave this topic here. 

1.3 A brief survey of the main results 

The results in this paper provide the foundations for the above Riemannian 
geometric approach to the three-body problem with zero angular momentum. 
In Section 2 we present the kinematic geometric framework for the reduction 
method which we shall work out, consisting of the two successive reductions 
M ^ M M* , where M ~ is the (congruence) moduli space and M* ~ S"^ 
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is the shape space. The second reduction uses the cone structure of M over M* 
to ehminate the scahng variable p = \fl by radial projection to the sphere. 

A trajectory in the configuration space M projects to its moduli curve 
7(t) in M, and away from the base point O the curve further projects to the 
shape curve 7*(f) on the sphere. Conversely, whereas 7(i) is determined up to 
congruence by 7(f), the real power of the above reduction method rests upon 
the knowledge of the subtle relationship which, in fact, generally exists between 
7(i) and the geometric (i.e. unparametrized) curve 7*. 

In reality, the complete study of the three-body trajectories is hereby re- 
duced to the study of the relative geometry between the shape curve 7* and the 
gradient flow of U* , namely the Newtonian potential function restricted to the 
2-sphere M*. The major results of the paper are divided into the following four 
main topics : 

• The calculation of the reduced Newton's equations in M in several ways, 
such as the geodesic equations of the Riemannian space {Mh,dsh), and 
the reformulation of the geodesic condition in terms of the curvature of 
the spherical shape curve. A suitable combination of these equations also 
yields a separation of the radial variable p and hence a third order ODE 
on the 2-sphere which describes all shape curves for any energy level h. 
The key step in this reduction is Lemma 3.8 which relates the intrinsic 
geometry of (7*, U*) to a kinematic quantity of 7(f). 

• The unique parametrization theorem (cf. Theorem 4.6) asserts that the 
time parametrized moduli curve 7(i) is (essentially) determined by the 
oriented geometric shape curve 7*. Furthermore, the curve 7* is in fact 
uniquely determined by the first two curvature coefficients at a generic 
point on the curve. The basic technique used here is the local analysis of 
solution curves via power series expansion. 

• The monotonicity theorem (cf. Theorem 5.8), which describes a type of 

picccwise monotonic behavior of the shape curve 7*. Namely, the mass- 
modified latitude is a strictly monotonic function along 7* between two 
succeeding local maxima or minima, and they lie on opposite hemispheres. 
In particular, the curve intersects the eclipse circle at a unique point be- 
tween two such local extrema. 

• Some initial applications to the study of triple collisions. In particular, 
simple geometric proofs of the fundamental theorems of Sundman and 
Siegel. Moreover, their asymptotic formulae for the derivatives of the 
moment of inertia / up to order two are extended to the derivatives of 
any order. 

The present exposition is based upon previous works of the authors (cf. [3], 
[4], [5]) on the three-body problem, exploring its "sphericality" as it manifests 

itself in various ways. The differential equations which describe the moduli 
curves of the three-body trajectories, are elaborated in Section 3 and 4, including 
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a careful power series analysis and comparison of the initial value problems at 
the moduli space and the shape space level. 

Section 5 is devoted to a geometric study of the gradient field of U*, which 
also yields a simple and purely geometric proof of the monotonicity theorem. 
A similar type of monotonicity for the shape curves was first proved by Mont- 
gomery [9] with his "infinitely many syzygies" theorem, and later by Fujiwara 
et al [2]. 

In Section 6 we recall the classical results and clarify some issues on the 
work of Sundman and Siegel concerning triple collisions. Moreover, with the 
results obtained so far, many challenging problems, for example in the study 
of collisions and periodic motions, naturally present themselves for an in depth 
study of the global geometry of shape curves. Some of these open problems will 
be briefly discussed in Section 7. 

2 Kinematic geometry of m-triangles 

A three-body motion with vanishing angular momentum is always confined to 
a fixed plane (for purely kinematic reasons), so the motions we shall study are 
always planary. Therefore, we choose a plane c with normal vectors 
±n and define an m-triangle to be a triple S = (ai, 32,33) of vectors a; G 
constrained by the center of mass condition in (8). Hence, for our purpose we 
shall modify the definition (8) of the configuration space by taking the subspace 

3 

M ~ C R*^ : ^ miBi = (16) 

i=l 

which consists of the above m-triangles in the fixed plane . M has the natural 
action of the rotation group S0{2), and the moduli space, representing the 
S'0(2)-congruence classes of m-triangles, is the orbit space 

M = M/SO{2) ~ RV5'0(2) d R^/S0{3) (17) 

The degenerate (or collinear) m-triangles constitute the eclipse subvariety 
E of M, and we say a nondegenerate m-triangle is positively (resp. negatively) 
oriented if (ai,a2,n) is a positive (resp. negative) frame of R^. Accordingly, we 
may write M and M as the disjoint union of three subsets 

M = M+UEUM- , M = M+UEUM_ (18) 

and moreover, we observe that the moduli space of unoriented m-triangles would 
be _ 

M/Z2 ~ IR''/0(2) ~ M± U £ ~ M.^/S0{3) = R^/0{3) 

Remark 2.1 For a study of general (non-planary) 3-hody motions, the natural 
configuration space M consists of pairs (6,rv), where 5 is an m-triangle in 
and n is a unit vector perpendicular to all a,. M is a 6-dimensional manifold (a 
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4-plane bundle over S'^) with the natural action of SO{3), and now the moduli 
space coincides with the above one (cf.[5], Section 2), namely 

M/SO{3) ~ M/SO{2) = M 

We shall describe in more detail the topology and induced Riemannian struc- 
ture of the above simple orbit spaces. Let 5 (as above) and 5' = (bi,b2,b3) be 
m-triangles. The following 50(2)-invariant, but mass dependent inner product 

(<S,<5') -^m,ai-bi (19) 

is just the kinematic metric ds^ of M defined by (10). In particular, the squared 
norm is the moment of inertia, = I = p^. Let 5^ C M be the unit sphere 
(p = 1) and denote its spherical metric by du^. Then we can express M as the 
Riemannian cone over {S^, dv^) 

M = ]R4 = C(S'3) : ds^ = dp^ + p^du^ (20) 

where p measures the distance from the base point (origin) of the cone. 
A description similar to (20) applies to the moduli space 

(M, ds^) = (M, ds^)/SO{2) (21) 

whose "unit sphere" M* = (p = 1) is called the shape space, namely it is the 
image of in M 

{M*, da^) = (S^ du'^)/S0{2) (22) 

with the induced metric denoted by da"^. Thus, M also inherits the structure 

of a Riemannian cone over its "unit sphere" 

M = C{S^)/SO{2) = C{S^/SO{2)) = C(M*) : ds^ = dp'^ + p^da^ (23) 

with p still measuring the distance from the base point O. The shape space is 
actually isometric to the 2-sphere of radius 1/2, as follows from the well-known 
Hopf fibration construction 

53 ^ M* = S^/SO{2) = CPi ~ 5^(1/2) (24) 

As a cone over 5^, M is clearly homeomorphic to with the origin at 
the base point O. Away from O they arc even diffcomorphic, when M has the 
induced smooth functional structure as an orbit space of M**. 

For the convenience of applying vector algebra we also recall the Euclidean 
model of M, where M is identified with M^, with Euclidean coordinates {x, y, z) 
and associated spherical coordinates {I,(f,9), and the kinematic metric is ex- 
pressed as the following conformal modification of the standard Euclidean metric 

^ dx^ + dy^ + dz^ ^ 2 ^ 2 ^ ^^2 ^f)2^ (25) 

= p'i = +y^ + z^, < <p < TT, < 6* < 277 
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Here (tp, 9) denotes any choice of spherical polar coordinates on the sphere 

whose induced metric from the Euclidean 3-space is that of the round sphere of 
radius 1 

5^(1) : = _^ gj^2 ^ ^^2 ^26) 
whereas its induced (i.e. kinematic) metric as a submanifold of {M,ds^) is 

da'^ = lrfs2 (27) 
By (15) and (25) the total kinetic energy can be written as 

T = f + T'^ = ]--p' + ^{^'^ + (sin^ ^ye") + T'^ (28) 
z o 

where the rotational term T'^ vanishes precisely when i7 = 0. Starting from 
Section 3 this is our standing assumption. 

Remark 2.2 The radial distance function p = \fl is also referred to as the 
hyper-radius in the physics literature. For our purpose it is generally more con- 
venient to use p rather than I as the scaling parameter, and we shall refer to 
{p, (fi, 9) as a spherical coordinate system of M. We refer to [5] for the relation- 
ships between spherical coordinates, individual moments of inertia {Ii, 12,13), 
or mutual distances (ri2, r23, rsi). 

In the above Euclidean model of kinematic geometry the decomposition in 
(18) has a distinguished equator plane, namely the eclipse plane E which divides 

into the two half-spaces M± = Mj_. We choose the Euclidean coordinates 
so that E is the xy-planc and the half-space z > represents the congruence 
classes 6 of the positively oriented m-triangles. Finally, 

E* =EC\M* -.x^ = l,z = Q 

is the distinguished equator or eclipse circle of the sphere M* = S^. 

On the other hand, the position of the various shapes (5*of m-triangles on 
the sphere is uniquely determined by the position of the three binary collision 
points hi, i = 1,2, 3, along the circle E*, where bi represents the shape of the 
degenerate m-triangle with a2 = as etc. We are still free to choose the cyclic 
ordering bi — > b2 ^ ba either in the eastward or westward direction. Moreover, 
the mass distribution and the relative positions of the three points b^ mutually 
determine each other. In fact, the angle between b2 and ba is given by 

mamg-mi 

cosp, = etc. (29) 

TO2m3 + mi 

and these formulae can be inverted. 
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Let {iV, 5} be the north and south pole of S"^. It is often convenient to 
choose the spherical coordinate system {<f, 9) centered at A^, namely if = Q 
at the pole A''. Then the poles represent the m-triangles (congruent, but with 
opposite orientation) of maximal area for a fixed size p, and more generally, the 
area of an m-triangle is given by the formula 

2 

A= . , ^ z|cos</p| (30) 

For a normalized (i.e. p = 1) m-trianglc of shape p e 5'^(1) we also recall 
the formula for the mutual distances (cf. (2)) 




P-bfel = J - — ^sino-fe (31) 
y niimj 

where |p — bj,| (rcsp. 2ak) is the Euclidean distance (resp. angle) between p 
and bfc. For proofs of (29) - (31) we refer to [4] or [5]. 

Remark 2.3 In this paper we use both the kinematic and Euclidean model 
S^{r), r — 1/2 or 1, of the shape space M* . Their arc-length parameters are 
a and s = 2a, respectively, cf. (27). Of course, the various geom,etric quan- 
tities, such as velocity, geodesic curvature, gradient etc. must also be scaled 
appropriately when passing from one model to the other. 



3 Analysis on the moduli curves of 3-body mo- 
tions with zero angular momentum 

In this chapter we shall follow Jacobi's geometrization idea at the level of the 
moduli space M = R^. This enables us to reduce the analysis of 3-body trajec- 
tories with zero angular momentum to that of the corresponding moduli curves. 
According to Jacobi, for a given energy level h the moduli curves can be in- 
terpreted as the geodesies of a specific Riemannian metric dsf^ on M. Now, 
the standard procedure for the calculation of the geodesic equations amounts 
to the calculation of the Christoffel symbols of the metric, with respect to a 
suitably chosen coordinate system suggested by the geometry of the space, say. 
The resulting equations are ordinary differential equations whose solutions are 
curves parametrized by the arc-length. 

The kinematic geometry describes M with the scaling and rotational symme- 
try of a Riemannian cone over a sphere, and therefore the spherical coordinates 
p,ip,0 present themselves as the most natural choice. In Section 3.1 we shall 
calculate the associated differential equations. However, the natural parameter 
for 3-body trajectories is the physical time t, and it is the effective usage of fixed 
energy that enables us to express the equations in terms of t as well, cf. (34). 

On the other hand, in the Hamiltonian least action principle the time interval 
is fixed and t is the natural parameter from the outset. The Euler-Lagrange 
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equations for the Lagrange function L on M arc calculated in Section 3.2, and 
this approach, in fact, yields the same system (34) in a much simpler way. 

However, we also seek a differential equation purely at the shape space level, 
that is, with the scaling parameter p eliminated. This demands a deeper im- 
derstanding of the relationship between p and the geodesic curvature of the 
shape curve. To this end we shall introduce an alternative geometric approach 
to the study of geodesies in M , which takes the full advantage of the spherical 
symmetry and the cone structure of M. This is the topic of Section 3.3. 

Finally, in Section 3.4 we shall synthesize the results obtained in the pre- 
vious subsections and, in particular, we explain how the moduli curve can be 
reconstructed from the shape curve. 

3.1 Calculation of the standard geodesic equations in 

We shall calculate the (standard) geodesic equations of (M^, ds\) relative to the 
spherical coordinates (p,(p,6), where {(p,0) denotes a choice of spherical polar 
coordinates on the shape space M* = S'^. The homogeneity of the Newtonian 
potential function (1) allows us to write 

U=-U*{^,0) (32) 
P 

where U* is the shape potential function on 5^. By (13) and (25), the metric 
with the arc-length element Sh is 

dsl = {— + h)ds' ^{— + h)( dp' + ^{d^' + sin2 ^ de')) (33) 
P P V 4 J 

and the standard procedure for the calculation of the geodesic equations via the 
Christoffel symbols, but expressed with respect to time t as the independent 
variable, yields the following system (cf. [13], (5.11)) 

(i) = p+p1-1{1u* +2h) 
P P P 

{ii) = ^ + 2^^ - ^ sin(2<p)^' - ^U; (34) 

(Hi) = + 2(^9 + 2 cot{(p)^e - ^^^m 

p p"* sin (f 

where equation (i) is just the Lagrange- Jacobi equation (5). 

The calculation of the above system goes as follows. For simplicity, let us 
write 

f = U + h = T, u = Sh, 
and d^/du = d^/dt = ^ for any function ^. Then by (15) and (33) 

du' = Tdt = 2T^dt'' = 2fdt^ 
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and we deduce the useful identities relating arc-length and time derivatives 

U U'^ U J 



^' = i, r = ^, -=24 (35) 



2f =p' + ^{'P' + (sin^ ^)^'), 2f = (36) 

Now, the first step is to calculate the geodesic equations with u as the in- 
dependent variable, following the standard procedure, and the next step is to 
transform the equations by changing to t as the independent variable, using 
the identities in (35) and (36). With the notation {xi,X2,X3) = {p,ip,9) the 
equations in the first step can be stated as 

= 0, * = I,2,3 (37, 



where the Christoffel sysmbols Vf^ are defined by 



dxi dxj dxr, 



Here (g*^) is the inverse of the matrix {gij) representing the metric du^, in the 
sense that 

du^ = P{dp^ + —{dip^ + siv? dO^) = ^^gijdxidxj 
Thus the matrices are diagonal and 

p2 ^ p2 2 /*2 2 • 2 77 —1 

511 = / ,922 = ^f 2:1,333 = ^/ a;ism X2, 5 = 9ii 

The first geodesic equation follows from (37) with A: = 1, so by calculating 
the symbols Vjj we obtain 

" ^* >P" + (S - P)%" + (sin^ 'P)0") + 4zW + U*eO') = 



Then, in the second step a straightforward calculation using (35) and (36) yields 

P+ ^ \ — TT^ - 2 ^ — + (72 2-)/3 ^ = 

P Pp^ P P^ PP f 



where the bracket expression simplifies to —p~^{U + 2h). The final result is the 
first equation of (34). 

Next, let us consider the case k = 2 where the first step yields the equation 

= ^" - 1^ {p'" + t(^" + («i^' ^)^")) + tI^" - ^(«i^2^)^" 

/ p p f p 
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The second step leads to the equation 



p 2 p'^ ^ 

f ■ ^u; -2 u* . u; 
/ f^p f^p^ f^p 



where the bracket expression simplifies to 

<p ,U* U* . d ^ 
^ -p--f) = 

P p r dt 

This yields the second equation of (34). The last case A; = 3 is similar to the 
previous one, so we have omitted the calculations. 

3.2 An alternative derivation of the standard geodesic 
equations 

Another way of deriving the system (34) is to calculate the Euler-Lagrange 
equations for Hamilton's least action principle (7), at the level of the moduli 
space M. In fact, the Lagrange function L = T + U descends to a function 
defined on the tangent bundle TM since n = 0, and 



1 

^=2 



d'y 
'dt 



2 



1 ds 2 

2^dt' 



is actually the kinetic energy of a curve ^{t) in M. Thus we have a simple clas- 
sical conservative mechanical system with potential energy —U, kinetic energy 
T, and conserved total energy h = T—U. Therefore, in terms of the coordinates 
(p, (f, 6) the associated Lagrange system is 



d dL _dL d dL _dL d dL _ dL 
'dt~dp~~dp' did^~d^' di'd§ ~ W 



(38) 



where by (25) 



L = T + U = l^p^ + ^{(p" + sin2 ^ + -U*{^, 



X ^ J P 

Now, the first equation of (38) reads 



dt'' 4''' ' " ' p 



-p=ti{^^+sm^^e)-^U* (39) 



and by substituting the expression 

^(V?2 + sin^ ^ 0^) = ^{2T - p2) = ^(2U + 2h- p^) 
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into (39) the equation transforms to equation (i) of (34). 
Next, the second equation of (38) reads 

and after differentiation this becomes 

p'^lp + 2ppCp - ^ sm2ip 0^ - -U* = 
I p 

which is precisely equation (ii) of (34). Similarly, the third equation of (38) 
reads 

-^(p2sinV6') = -f/; (40) 
at p 

where the left side becomes 

Ipp svT? (f 6 + 2p^ sin ip cos ipO + p^ sin^ (p 
Then it is easily seen that equation (40) becomes equation (iii) of (34). 

3.3 Cone surfaces and geodesies in 

To take full advantage of the cone structure of M and the scaling property (32) 
of the potential function U, one naturally seeks to reduce the analysis of the 
moduli curve to that of the shape curve. Therefore, in this subsection we shall 
study the geodesic equations with the shape curve in the forefront. 

Definition 3.1 Let j be a curve in M not including the (triple collision) base 
point O. The cone surface spanned by 7 consists of all rays emanating from O 
and intersecting 7. The cone surface is denoted by C(7) or C{j*). 

The intersection of the cone surface with the shape space M* is the associ- 
ated shape curve 7* , and conversely, the cone surface is also uniquely determined 
by 7*, which explains the notation C(7*). 

Let 7(t) be the moduli curve of a 3-body motion with zero angular mo- 
mentum and 7*(i) the associated shape curve with arc-length parameter a in 
M* ~ 5^(1/2). Endowed with the kinematic metric, the cone surface 

C(7*) : ds'^ = dp^ + p^da'^, < < ai (41) 

is isometric to a (flat) Euclidean sector with polar coordinates (p, a) and angular 
width a I. 

Next, let the function u{a) denote the restriction of U* along 7*. There is a 
pair of conformally related metrics on the above sector, namely the flat metric 
(41) and the dynamical metric 

dsl = {h + ^)ds^ (42) 
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defined on the subregion 



Ch*)h = {{p,'j);p>0,0<a<ai,h+^>0} (43) 



On the one hand, 7 is a geodesic curve of the surface {C{"f*)h with the 
metric (42), but on the other hand, it is also a geodesic of the ambient space 
{Mh, dsf^) and hence the normal component of the curvature vector of 7 in 
also vanishes. Accordingly, it is natural to write the geodesic equations of 7 in 
Mh as a pair of coupled ODE's expressing, respectively, the vanishing of the 
tangential and normal component of the curvature vector. 

To analyze curvatures, let us fix some convention concerning orientation. 
We assume 7 (and hence also 7*) is oriented; they are curves in M ~ and 
M* ~ S"^ respectively, and these spaces have their standard orientation. We 
choose a positive orthonormal frame (r,r7,i/) of (Af, ds^) along 7, as follows. 
Let T (resp. r*) be the positive unit tangent field of 7 (resp. 7*), and choose 
the unit normal field f* of 7* so that (r*, i^*) is a positive frame of M*. Then 
1/ = {1/ p)h>* is a unit normal field of the cone surface (7(7) and hence orients 
the surface. Finally, r) is the normal field of 7 in C(7). For convenience, we 
write 

d . I d . d Id 

T = cosa— — hsma-— — , 7? = — sma— — hcosa-— (44) 
op poa op p Oa 

where a denotes the angle between the radial and tangential direction, that is, 
the angle between ^ and t. 

Let n be a vector normal to 7, say rj 01 v as above, and let /C(n) and /C?i(n) 
denote the associated geodesic curvatures of 7 with respect to the metrics ds^ 
and dsj^, respectively. It follows from the first variation formula of arc-length 
that the curvatures of a given curve with respect to such a pair of conformally 
related metrics are linked by the following formula 

}Ch{n)=}C{n)~-^ln{h + U) (45) 

There are two cases to analyze, namely n = 77 or i^, and for simplicity we write 

/C — IC{ri), K.^ = ^{^) etc., and we assume 7* is not a single point. 

First, the geodesic curvatures of 7 in the surface C(7*) with the two metrics 
are related by 

^ d , ,da da. Id,,, u(a) . 

/C. = /C - - - ln(. + U) = (- + -)--- Hh + ^) (46) 

where we recognize (a + a) as the angle between the tangent line and a fixed 

reference ray in the Euclidean sector (41). Similarly, let /C^ and /C^ be the 
"surface normal" geodesic curvatures of 7 with respect to the metrics ds"^ and 
dsf^, consequently 
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To find an expression for /C^, it is a key observation that the principal 
curvatures of C(7*) in M at a given point {p,a) are the numbers and jfC* , 
where /C* = JC*{a) denotes the geodesic curvature of 7* on the sphere M* ~ 
5*^(1/2) at the point 7*(o"). Therefore, by the classical Euler's formulas, the 
normal sectional curvature of C(7*) in the tangential direction of 7 is equal to 

IC^ = ^{sm'^a)IC*{a) (48) 

Finally, the geodesic condition for 7 in reads 

which by (46)-(48) is neatly expressed by two scalar ODEs. In terms of the 
kinematic arc- length parameter ds of 7, regarded either as a curve in C(7*) or 
M, we can state the final result as the following theorem : 

Theorem 3.2 Let ^{t) (resp. j*{t)) be the moduli (resp. shape) curve of a 
given three-body motion with zero angular momentum and total energy h. Set 
a to he the arc-length parameter of in the shape space Ad* ~ S'^(l/2), (p, cr) 
the polar coordinate system of the associated cone surface C(7*), and u{a) the 
restriction of U* along 7* . 7/7 is not a ray solution, then it is characterized by 
the following pair of equations 

da da 1, . d cos a 9 , , u{(t) ^ ^ 

{ii): {sm'a)]C*-l^Hh+^)=0 (50) 

where a is the angle between the radial direction and the tangential direction 
and V* is the positive unit normal vector field of 7* in M* . 

Remark 3.3 The exceptional case of ray solutions, that is, 7* is a single point, 
can be settled directly from (45), where /C(n) = /C/i(n) = and hence for each n 

\n{h + U) = or Vt/ • n = 

an 

This is equivalent to VC/* = 0, so 7* is one of the five critical points of U* 
on the 2-sphere, namely the two minima (called Lagrange points) and the three 
saddle points (called Euler points) lying on the equator (or eclipse) circle. 

We also deduce the above result from the equations (ii), (Hi) of (34), namely 
U* and U0 must vanish if (p and are set to be constant. Ray solutions yield 
the simplest type of three-body motions, namely the shape invariant or so-called 
homographic motions, which in the case of ft = are confined to a line. Then 
the only variable p{t) is the solution of the 1-dimensional Kepler problem given 
by the Lagrange- J acobi equation (i.e. (i) of (34)). 
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Remark 3.4 (i) For a givnn value of total energy h, the influence of U* on 
the geometry of the associated cone surface C{j*)h is via the function u = 
?7*|-y. , and equation (49) is exactly the geodesic equation of C{'y*)h- Most of 
the geodesic curves ofC{'j*)h are, of course, not geodesies of the ambient space 
Mh since they are not moduli curves of actual three-body motions. 

(ii) Equation (50), on the other hand, is expressed in terms of the relative 
geometry of the inclusion 7* C M* , namely the geodesic curvature K* and the 
normal derivative of U* along 7* . 

Finally, recall from (44) that tlw^ sc;aling variable p and the angular variable 
a, which measures the radial inclination of the moduli curve in the cone M, 
essentially determine each other via the relations 



which tell us, for example, how to calculate p from the variation of a along the 
shape curve 7* C 5^(1/2) : 



3.4 Synthesis of the analysis of the moduh curve and that 
of the shape curve. 

In the previous sections we have used geometric ideas to obtain differential 
equations in the moduli space M ~ characterizing 3-body trajectories with 
zero angular momentum. Since M is a cone over the 2-sphere M* defined by 
p = 1, it is natural to project the moduli curve 7 down to its image curve 7* 
on the sphere. However, unless one resolves the hidden interlocking between 
7* and the scaling variable p, implicitly described by the differential equations, 
one cannot reconstruct the moduli curve from its shape curve and thus fully 
utilizing the reduction from M to M* . 

In order to separate the scaling variable p from the spherical variables (tp, 9) 
we shall proceed by combining the two systems of geodesic equations, (34) and 
(49) - (50), which we derived in two different ways. Note that the natural 
parameter in mechanics is the time t, whereas the arc-length parameter is the 
natural parameter in metric geometry. This suggests a transformation of the 
latter system to equations with t as the independent variable, and as it turns 
out, this also provides a remarkable simple solution of the above separation 
problem. 

3.4.1 Basic geometry of curves on the 2-sphere with a potential func- 



Since spherical curves play a crucial role in the present study, it is convenient 
to collect some basic formulae concerning the differential geometry of curves on 
5^(1), as well as the tangential and normal derivatives of a given (potential) 




(51) 



P((t) = p{aQ)e^''o <:o*«(<^)<^'^^ for p(a-o) ^ Q 



(52) 



tion 
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function U* on the sphere. We shall express them in terms of a chosen spherical 
polar coordinate system {(p, 9). 

For a given oriented curve 7*, let r* (resp. u*) be the unit tangent vector in 
the positive direction (resp. unit normal vector) such that (t*, u*) is a positively 
oriented frame of the sphere. We consider a (regular) time parametrized curve 
7*(t) = {cp{t),9{t)) and set s = s{t) > to be the arc-length along the curve. 
As before, differentiation of a function / with respect to t or s are denoted by / 
and /' respectively, and clearly / = f'v where i; = s is the speed of the curve. 
Then 

r* = ^ = i(^f +4), 1^* = V^sin^A+^^|) (53) 
ds V dip oO V dip sm ip o9 

and the velocity vector field of the curve is 

^=vr*, v = \[^^^(^^ (54) 

The scalar acceleration 

V = —V = —[ifip + {sinp cos ip)ip9 + sin^(<^)^^] (55) 

and its higher time derivatives are needed to express time derivatives of a func- 
tion in terms of arc- length derivatives, using operators of increasing order 

dt=''d-s' dt^='d^-'''d-s^'''- ^'^^ 

To calculate the geodesic curvature function K*, let us first make use of 
Euclidean coordinates 

X = sinp cos 9, y = sin ip sin 9, z = cos (p 

and write x(s) = {x{s), y{s), z{s)) and use the formula 

i<:*(s)=x(s) xx'(s)-x"(s) 

where x'(s) = r*. This yields 

K* = {cosip)9'{l + p'^) + smp{p'9" - 9'p") (57) 
= ^ ^{cosp)9{v'^ + ip^) + s\up{>p9 - 9!p)^ 

and its intrinsic first derivative is 

K*' =-^K* = {-{sinp)p'9' + {cosip)0"){l + ip'^) + {cosp)ip'{ip'9" - 9'p") 
+ sinp{p'9"' -9'p'") (58) 
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The gradient field of U* is the following vector field on the sphere 



d_ u; d 



WU* = Ul— + -^- (59) 



which allows us to calculate various derivatives of U* . For example, the tangen- 
tial and normal derivatives along the curve are, respectively, 

dU* 1 

= ^ = ^u* • T* = - {cpu; + eu;) (eo) 
c^: = = vc/* • = -{-esin^u; + (p-^u;) (ei) 

ou* V ^ smip 

and the intrinsic first derivative of U* is 

U:' = = i-^W + {-sin^U;)e" + (62) 

as sm (fi ^ sin ip 

+ (4^ - u;e sm^)e'^ + {-^ - sm<pu;^ - cosipu;)<p'e' 

simp ^ sijup ^ 

Finally, for convenience and later reference let us introduce the following 
definition : 

Definition 3.5 For a given curve 7* and function U* on S^, the associated 
Siegel function along the curve is defined to be 

U* 

S = # (63) 

We regard the function as undefined along geodesic arcs. Moreover, the 
function may have a singularity at isolated points where K* vanishes. Note 
that S is independent of the orientation of 7*. Observe the following formula 
for the logarithmic derivative of S, as a function of s (and similarly for t as 
parameter) 

6' U*' K*' 

^ = ^-fl_ (64) 

Remark 3.6 We have named the above function after C.L. Siegel for the follow- 
ing reason. In his study (cf. [11]) of triple collisions in the three-body problem, 
Siegel investigated the asymptotic behavior of the time derivatives I , I of the 
moment of inertia I = . The major step in his proof was, indeed, to show that 
the expression VT(2r - P/4:T) tends to zero. It turns out that this expression 
equals & whenever the latter is defined (see Lemma 3.8, where (f'v^ equals the 
above expression) . In particular, it is an intrinsic quantity at the shape space 
level. We shall return to triple collisions and Siegel's approach in Section 6. 
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3.4.2 Reformulation of the geodesic equations in terms of the shape 

curve 

For convenience, let us write = ^ and V2 = ^ for the basic coordinate 
vector fields of the spherical coordinate system on the unit sphere. By definition, 
the acceleration of the spherical curve 7*(f) = {ifi{t),9{t)) is the expression 

which is the covariant derivative of the velocity along the curve, with respect to 
the metric (26). By definition, 

and the only nonzero Christoffel symbols of the metric are 2 = = cot f, 
T 22 = — sin If cos (p. Consequently, 

• 2 d •• • d 

^* = (m — 9 sin w cos (p) 7- — \- (0 + 2(h0cotip)— (65) 
oif oO 

Now, take the above sphere and curve to be the shape space M* and a shape 
curve 7*(t), respectively. Then it follows immediately from (59) and (65) that 
the equations (ii) and (iii) of the system (34) can be expressed neatly as the 
following coordinate-free vector equation on the 2-sphere, 

7* + (2^)7* - ^VC/* = (66) 

Here the scaling variable p = \fl of the cone M plays the role of an auxiliary 
function which couples equation (66) to equation (i) of (34). The latter is 
the Lagrange- Jacobi equation (5), which for a given shape curve 7*(t) is a 
second order differential equation purely for p{t). Another interpretation of the 
coefficient of the velocity in (66) follows from (51), namely we have 

2^=f;cota (67) 
P 

Remark 3.7 In j act, the differential equation (66) can he completely decoupled 
from the function p{t) and hence it really becomes a differential equation on the 
2-sphere, as summarized in Theorem 3.9 below. 



3.4.3 Separation of the scaling function p{t) from the shape space 

coordinates 

The system (49)-(50) characterizes the moduli curves of 3-body trajectories with 
zero angular momentum, expressed in the language of kinematic geometry and, 
in particular, for that reason the natural parameter is the arc-length s of the 
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moduli curve or the arc- length cr of its image shape curve on the sphere of radius 
1/2. Although the original mechanical system (2) naturally involves the physical 
parameter of time t, the latter is infinitesimally related to s by the identity 



(68) 



and this enables us to express the system (49)- (50) in terms of t and hence 
combine it directly with the other system (34). 

To this end, let us consider a shape curve 7*(f) on the sphere 5^(1/2). By 
(51), 

, = [ n. H I " ' ~ 

dt' dt 

and then equation (50) becomes 

2 

IC* = —-—U* 
2p du* 



da ds p t^*N-i/2 
sma = p— /— = — =(/i H ) ' 



dt 



El 

2 



d7* 



dt 



1 d 



(69) 



(70) 



The latter is not only considerably simpler than equation (50), but it also pro- 
vides a simple formula to compute p{t) — p{'^{t)) in terms of the geometry of 
the shape curve, namely 

(71) 



IC* 



dt 



In view of the integral formula (52) this is, indeed, a pleasant surprise which, 
in one stroke, shows how to reconstruct the moduli curve ^{t) from the shape 
curve 7*(t) by the simple formula (71). 

The expression on the right hand side of (71) refers to the kinematic geometry 
with M* = 5^(1/2), and the whole product on this side would change by the 
factor 4r^ if we had worked in the sphere S"^ (r) . Henceforth, we shall return to 
the unit sphere <S'^(1), and with the notation for speed, curvature and normal 
derivative from Section 3.4.1 we can restate (71) in the following way : 

Lemma 3.8 Let & be the Siegel function (63) of {j*,U*), which relates the 
intrinsic geom,etry of 7* with the gradient field VU* on the unit 2-sphere. If 
7(t) = {p{t),"f*{t)) is the time-parametrized moduli curve of a three-body motion 
with zero angular momentum, then the speed v{i) of j*{t) is related to p{t) and 
6(7*(t)) by the identity 



or 



S = ^p^v-" 



(72) 



In particular, & is always nonnegative! 
Now, returning to equation (66) we set 



_ 2 6 4w 
p~3&~3v 



Q = - 



V 



(73) 
(74) 
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where the rightmost identity in (73) foUows from (logarithmic) differentiation 
of the identity (72). Then we can state the foUowing result : 



Theorem 3.9 For 3-body motions with zero angular momentum and fixed total 
energy, the associated shape curves 7* (<) on the unit 2-sphere are characterized 
by the ODE 

^,*+P^p*+QVU^=0 (75) 

where the first term is the (covariant) acceleration and the coefficients are the 
functions P,Q defined by (73), (74), which can be expressed purely in terms of 
tp, 9 and their derivatives up to order 3. 

Remark 3.10 The formula (72), or equivalently 

4U* 

K* = ^ (76) 

expresses the geodesic curvature of the shape curve in terms of p,(p,9,(p,9. 
Therefore, since it only involves first order derivatives in the moduli space, it 
is not surprising to find that the same formula can, indeed, be derived more 
directly from the general spherical curvature formula ( 57) by elimination of the 
second order derivatives using equations (ii), (Hi) in (34)- 

3.4.4 Regular and irregular points and exceptional shape curves 

The formula for in (72) involves the three quantities v, U* and K* and the 
product may become indefinite when some of them vanish, namely we note the 
following implications 

[v^0]=> [U* 0] ^ [K* = 0] (77) 
[C/; = 0] =^ [v = 0] or [K* = 0] 

It is worthwhile having a closer look at the geometric interpretation and behavior 
of the shape curve due to the vanishing of any of these numbers, and accordingly 
we shall make some definitions to distinguish the various cases. 

Formula (72) expresses the Siegel function 6 of 7* in two different ways, 
namely at a point Pq = 7*(io) its value is 

limjA. = eo=Hm§ (78) 

whenever any of the limits are defined, including +00 as a limiting value. 

Definition 3.11 Pq is a regular point if < &o < 00, and otherwise it is 
irregular. 
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The irregular points consist of the cusps, the collision points (binary or 
triple), and for completeness we also include escape points : 



cusp at Pq 
triple collision at Pq 
binary collision at Pq 



60-0, «o = 0, V;7*(Po) 7^0 
60-0, Po=0, Vf/*(Po) = 
60 = 00, vo = 00, Pq = bi 



(79) 



escape at Pq 



So =?, Po = oo, vo = 



In Section 6 we shall return to triple collisions, but binary collisions and escape 
to infinity behavior will not be a topic in this paper. 

The simplest curves on the 2-sphere are the geodesic circles, characterized 
by K* = at each point. If such a curve is the shape curve of a three-body 
motion, then it follows from (78) that U* also vanishes. In fact, 7* coincides 
with a gradient line segment if and only if it is a geodesic. 

Definition 3.12 The shape curve 7* is called exceptional if it is confined to a 
gradient line (or a geodesic circle), or it consists of a single point. 

Clearly, a single point shape curve must be a fixpoint of the gradient flow, see 
Remark 3.3. Thus, apart from the exceptional shape curves, being a regular or 
irregular point is an intrinsic property, that is, it depends only on the geometric 
curve. 

Remark 3.13 We omit the proof here, so we rather claim that the only excep- 
tional shape curves (of length >0) are those representing collinear motions or 
isosceles triangle motions. Their crucial property is that the shape curve lies on 
a circle fixed by an isometry (reflection) of the sphere which leaves U* invari- 
ant. It also follows that a non- exceptional shape curve 7* can only intersect an 
exceptional curve transversely, that is, neither tangentially nor with zero speed. 

The equator circle E* represents the collinear motions, of course, but isosce- 
les motions exist only for special mass distributions, as follows. An isosceles 
m-triangle has (at least) two equal masses, say rni = TO2, and the mass ms 
lies on the symmetry axis of the triangle. Their shapes constitute the meridian 
through the north pole, the Euler point es and its antipodal point — es = bs. 
The latter point represents the collision of the two symmetric mass points some- 
where on the symmetry axis. Thus, an isosceles triangle motion arises when the 
initial position and velocity have the above isosceles symmetry. 

Henceforth, we shall assume the shape curve 7* is not of exceptional type, 
unless otherwise stated. Consider the power series expansions 



at a regular point Pq ~ 7*(to), where s is the arc- length measured from Pq. The 
value So of 6 at s = can be calculated in two ways, possibly by the aid of 



00 



00 
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I'Hospitals rule, 



Ion U* lim-fU* lim^rC/* Uk 

-pIv^ = 6o = lim = = ... = = ^ > 80 

where A; > is the smallest integer such that K/. ^ 0. Then Ki = Ui = for 

i < k, LUk ^ 0, and wc say Pq is a regular point of order /c. It is a finite number 
since otherwise U* and K* vanishes identically and 7* would be exceptional. In 
particular, we see that 7^ means 7* (t) is transversal to the gradient flow at 
t = to. 

The order fc of a cusp at Pq = {(p^, 9q) can be defined similarly by considering 
the limits in (80). The only difference is that Wfe = and ^ at the last 
step. Cusps arise when the moduli curve ^{t) in M is tangent to the ray at the 
point (pq,Po) and hence the projected curve 7*(t) on the 2-sphere "halts" at 
t = to- Geometrically, the curve 7* near Pq is a cusp consisting of two diverging 
branches which emanate from Pq, both with the initial direction of VJ7*(Po) 
and the initial curvature 

^o = i^ln|VC/*(Po)| (81) 

In Remark 3.17 we further describe the totality of cusps at Pq, at a fixed energy 

level h, as a specific family of curves 7* parametrized by a number c > (resp. 
c > 0) for ft > (resp. h < 0). In the case c = the two branches coincide 
completely and 7(^0) = {Pq,Po) lies on the Hill's surface dM^, cf. (14). 

3.4.5 On the initial value problem for the moduli curve and the 
shape curve 

In the spherical coordinates (p, if, 0) on M ~ M^, the time parametrized moduli 
curve and its associated shape curve are simply related by 

and conversely, formula (72) is the key to the lifting procedure, namely the 
reconstruction of "f{t) from its projection on the 2-sphere. 

We choose an initial point Pq — {fo^ ^0) on the shape curve, and assume (for 
simplicity) it is regular in the sense of Definition 3.11. In particular, the (local) 
lifting procedure is well defined, and the corresponding initial value problems 
for the ODE's (34) and (75) respectively, are equivalent. Thus, on the one hand, 
the solution ^{t) is uniquely determined by its initial position and velocity 

7(io) = (Po> "Po^ ^0), ^7(*o) = (Pi, Vi, Oi), (82) 
but on the other hand it is also determined by the corresponding initial data 

(¥'o>6'o),(¥'i,6'i),(<P2,6'2), where ip2 = ^ip\to, 02 = ^0\to, (83) 
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at the shape space level. 

However, the above chosen initial data determine an energy level h, and 
conversely, if h is already given, the initial data (82) or (83) are acceptable only 

if the resulting energy level is h. To make this more transparent, let us rather 
state the initial value problem in the moduli space as a system with 4 equations 





= 


p+ — --(-£/* +2/i) 
P P P 




{ii) 


= 


if + 2^(p - ^ sm{2(fi)e^ 
P 2 


4 

- —U* 


{Hi) 


= 


9 + 2-0 + 2cot(ip)ipe - 
P 


4 1 . 
p^sinV ' 


{iv) 


= 







where the first order equation (iv) is the energy integral (4) for the fixed value 
h. In fact, any one of the four equations is redundant and can be deduced from 
the other three. For example, with h calculated from the initial data (82), the 
solution of (i)-(iii) also satisfies (iv). 

More geometrically, if the initial data set (82) is a point on a specific level 
surface of type (iv) in the tangent bundle of M, then the solution j{t) of (i)-(iii) 
in (84) must lie on the surface (iv) for all t. But we can also determine the same 
solution 7(t) from (ii)-(iv). All this amounts to saying that for a given value of 
h the whole system (84) is of total order 5. 

On the other hand, the ODE (75) on the 2-sphere is independent of h and its 
total order is 6. Its solutions are the projections 7* of all solutions of the system 
(84) for any value of h. In general, they arc divided into three disjoint classes, 
distinguished by the sign of h (positive, zero or negative). The exceptions to 
this subdivision are precisely the exceptional shape curves (cf. Definition 3.12), 
which can represent three-body motions at any energy level. Sec Section 4.2. 

Now, let us consider the problem of how to translate the initial data (82), 
at a given energy level h, into a set of suitable initial data consisting of five 
numbers depending only on the shape curve. A natural first choice would be 

(<y?o,6'o),(V'o.f^o),So (85) 

where the angle f/'o ^ [0, 27r) specifies the initial direction of 7*, fo > is the 
inital speed, and So = {U^/K*)o > is the initial value of the Siegel function. 

Remark 3.14 The generic points of j* are the regular points of order k = 0, 

and hence 7* is transversal to the gradient flow almost everywhere. For such 
points it is, perhaps, more natural to replace So by Kq in (85). In fact, the two 
choices - either Kq or So - are equivalent since the direction tl)^ determines wq 
when k = 0. 
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The pairs (ipQ,vo) and ((pi,9i) evidently determine each other, namely they 
specify the initial velocity of 7*. Clearly, the 5-tuple (85) is merely the data 
(82) with the pair (pg, p-^) replaced by the single number &o, which is a second 
order quantity at the shape space level. Indeed, we recover immediately from 

(85) by making use of (80). 

Thus, it is clear that the initial data information given by (82) or (85) would 
be equivalent if we could also recover the radial speed p^^ from (85). However, 
the energy integral (iv) of (84) determines only pi; in fact, there is no solution 
at all if h is below the critical value 

hmin = IpIv'^ - —U*{(po, 0o), where p^ = {^^Y^^ (86) 
8 Pq 

There is a unique^ solution if h = /imin; with p-^ = 0, and for h > /imin there are 
two solutions which are distinguished by the sign of pi . 

A slightly different approach is to combine the identities (67) and (69), which 
yields a value of sin a, where < a < tt. In fact, a < 7r/2 means < and 
a > 7r/2 means Pi > 0. But, we are still left with the problem of how to 
determine the sign of Pi from the initial data (85). Anyhow, up to now we have 
the following result as a summary of the above local analysis at a regular point. 

Proposition 3.15 Consider the three-body motions with zero angular momen- 
tum and a given total energy h, whose oriented shape curve at a given regular 
point have the same initial direction, speed, and curvature (or Siegel number 
Sq; if the curvature Kq vanishes). Then the number of solutions is, up to con- 
gruence, equal to 0, 1 or 2 depending on whether h < /imin (resp. h = ft-min or 
h > /imin)> where /imin is calculated from the given initial data by the formula 

(86) . 

A three-body motion is said to be expanding when p(t) > 0, and it is con- 
tracting when p{t) < 0. Accordingly, we say the expansion index at time to is 
the sign e{to) = 0, ±1 of p{to) = Pi- 

Corollary 3.16 A three-body motion "/(t) with zero angular m,om,entum, and a 
given total energy h is uniquely determined up to congruence by the oriented 
shape curve 7* C S"^ as a subset (i.e. non-parametrized), together with the 
initial speed vq and expansion index e{to) at a regular point 7*(to)- 

Remark 3.17 The above corollary is, in fact, also true when the point Pq = 
7* (to) is a cusp. Indeed, for each fixed h, the family of possible cusps at Pq is 
parametrized by the nonnegative numbers 

c = p\ = 2{^^+h) (87) 
Po 

as follows from the energy integral (iv) in (84). 

For c = there is a unique moduli curve j{t) starting out from the "rest 
point" (poi-Po) on the Hill's boundary dMh, and the corresponding shape curve 
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7* is the simple cusp (i.e. with one branch) emanating from Pq. In general, the 
local geometry ofj* at Pq determines the number c. In fact, one can determine 
c from the first curvature coefficients Ki,i > 0, but note that Kq is independent 

of c, by (81). Finally, one solves the initial value problem in M at (pg, Pq), with 
the velocity component — ±\/c selected according to our choice of e{to). 

For a given non-parametrized curve 7* in M*, the geodesic curves of the 
associated cone surface C{'y*)h are characterized by the ODE of (49). The 
squared speed of such a curve is, of course, given by 2T = 2{U + h) and hence 
it is determined by the position. A geodesic 7 is therefore uniquely determined 
by the initial position and direction aXt = to, namely p{to), a{to) and the point 
7* (to)- By (71) or (72), and elimination of the speed in the shape space using 
(69), it follows that the validity of the identity 

. 2 26 

sin^ a = — — (88) 

along the entire curve 7, where 6 is the Siegel function of 7* in 5^(1), is a 

necessary and sufficient condition for a geodesic of the cone surface to be a 
geodesic curve of M/j as well. This proves the following statement : 

Corollary 3.18 A non-parametrized (i.e. geometric) curve 7* C <S''^(1) can 
be suitably parametrized as the shape curve of a three-body motion with zero 
angular momentum and total energy h if and only if its cone surface C{'y*)h 

has a geodesic curve saMsfying (88) along the entire curve. 

Problem 3.19 As indicated by the condition (88), only a very special kind of 
geometric curves on the 2-sphere can be suitably parametrized as the shape 

curve of a three-body motion as above. How can they be characterized in a neat 
way? Can such a curve have different time parametrizations as the shape curve 
of three-body motions? 

Problem 3.20 Is the initial speed and expansion index also determined by the 
geometric shape curve in Corollary 3.16? (The case h = turns out to be 
special.) 

4 On the analysis of moduli and shape curves 
via power series 

We continue to use the notation and terminology from Section 3. Consider 
a time parametrized moduli curve t ^{t) in M ~ which represents a 
three-body motion with vanishing angular momentum, and let < — > 7* {t) be the 
associated shape curve, namely its projection in the 2-sphere 5^. In this section 
we shall investigate the possibility of reconstructing the parametrized curve 7(f) 
solely from the oriented geometric shape curve. Moreover, there is the question 
of how much geometric information about the curve 7* is really needed for such 
a lifting procedure. At the end we shall also answer the question concerning the 
uniqueness of the time parametrization. 
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4.1 Generation of recursive relations and intrinsic geomet- 
ric invariants 

In the local analysis of the moduli and the shape curve, and their interaction 
with the potential function U*, we shall distinguish between two types of vari- 
ables or quantities. Namely, on the one hand there are the intrinsic quantities 
which depend only on 7* as an oriented geometric (i.e. unparametrized) curve 
and U* as a function on S^, and on the other hand there are the variable quan- 
tities, defined along the curve 7 or 7*, which depend on the scaling function p 
in the moduli space M or the time parametrization of the curves. The basic 
intrinsic quantities are the gradient field VU* (or the tangential and normal 
derivatives U*, U*), the unit tangent field of 7*, and the geodesic curvature 
function K* of 7*. Moreover, we shall assume that 7* is not exceptional and 
hence the linkage between 7* and U* is also neatly encoded into the intrinsic 
Siegel function 6 = U*/K*, see (63) and Section 3.4.4. 

Let s be the arc-length parameter of 7* measured in the positive direction 
from a given regular point Pq of order fc > 0. Then the coefiicients of the power 
series expansions of the above functions, such as 

K* =Ko + KlS + K2s'' + ... 
U* = uo+uis+U2s'^ + ... (89) 
U* =oJo+ Wis + ijJis^ + .... 
S = So + 6is-F62S^ + ... 

yield intrinsic quantities localized at the point Pq- Note the expansion of the 
tangential derivative of U* is 

V* = \/U* ■ T = ^U* = Ml + 2U2S + 3U3S^ + ... (90) 

as 

and the coefficients &n are expressible as rational functions of Ki and u)i,k < 
i < n + k. Let us say the order of a coefficient is the highest order of derivatives 
of local coordinates in its expression. Thus, the pair ((^q.^o) and uq arc the 
intrinsic (geometric) data of order at Pq- Next, the triple ui, lvq, and the unit 
tangent vector at Fq represent the intrinsic data of order 1 at Pq, and Un and 
Un+i (resp. Kn and 6„) has order n + 1 (resp. n + 2). 

We choose a spherical polar coordinate system {<f,0), with Pq difi'crcnt 
from any of the "poles" tp = or it, and for a given moduli curve 7(t) = 
{p{t), (p{t), 9{t)) we shall expand the coordinate functions, as well as U* and its 
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partial derivatives, as power series with respect to t : 

P = Po + Pit + P2*^ + Ps*^ + •••• 
f = ^0 + fit + ^2*"^ + fat^ + •••• 
g = gQ + eit + 02t^ + Ost^ + 

V = Vq + Vi + V2t + V2t^ + (91) 

[/* = MO + Ult + U2f + mt^ + .... 
U;=I^O+ jJi^t + Il2t^ + IJL^t^ + .... 

C^e* =%+??l^ + %^^ + %^^ + ••• 
For convenience, we also write 

sin(2^) = /o + /it + /2i^ + .... 
s\T?{ip) = 50 + 51* + 52*^ + •••• 
and list some of the initial coefRciens : 

uo = U*{ipQ,eQ), ui = /io</5i etc. 

/o = sin(2(/?o), /i = 2cos(2(^o)(^i, etc. (92) 
go = sin2(<Po), 9i = /oVi, etc. 

vi = — [2(p^ip2 + sin(yo) cos(<^o)'^i^i + 2sin^(v?o)6'i6'2] 

Vo 

where the expression for vi follows from (55). We shall regard Ho, rjo as intrinsic 

data, but they depend on the coordinate system, of course. 

Below we shall investigate dependence relations among the coefficients Pj, 
iPj, 6k of the coordinate functions in (91) and various other coefficients. Some of 
them are directly expressible in terms of the intrinsic data and hence regarded 

as constants, whereas the others are the variables. 

Definition 4.1 The following list of coefficients from (91) 

p„,vo;Pi,LPi,9i;p^,(f2,d2 (93) 

will be referred to as the variables of order < 2. The variables of order n are 
Pn,^n^ whcn u > 0, and Po, vo are the only variables of order zero. 

Henceforth, assume the above moduli curve j{t) is a solution of the ODE 
system (84). By inserting the power series into the equations (i)-(iv) of (84) and 
applying the method of undetermined coeflicients, we arrive at the following 
scheme of recursive relations for the variables of increasing order 0, 1, 2.. : 

£^10 : = 2pIp2 + Pop\ - 2/iPo - '"o 

E2o:Q = 2pl^2 + "^pIpiVi - \poh^\ - 4^0 (94) 
£30 : = 2gople2 + 2ffoPoPi^i + plhfi^i - 4% 
-E40 : = Pqp\ + ^plifl + 9oOl) - 2uo - 2hpo 
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£■11 : = 6po/93 + 8P0P1P2 + Pi- '2'hpi - wi 

£21 : = 6/9^^3 + 10/9^pi^2 + 4(PoP2 + PqpD^i - 2/oPo^i^2 

-^(/iPo + 3/oPoPi)^?-4Mi (95) 
-E31 : = 6ffoPo^3 + (IO50P0P1 + 251P0 + 2/oPoVi)^2 + 2/op^^iV'2+ 
(/iPo + 3/oPoPi)¥'i^i + (450P0P2 + 4fifoPoPi + 25iPoPi)^i - 4r?i 
and in general 

£i„:0 = (n + 2)(n + l)p^p„+2 + 

£;2„ : = (n + 2) (n + l)p^ v'„+2 + (96) 

£^3n : = (n + 2)(n + l)p^6i„+2 + 

where the remaining terms are of less order since they involve Pj, <^j, Oi for i < 
n + 2. For example, the coefficients u„, /i„, 77^ occur in (96) and their order is 
n. The equations £^4„ for n > are omitted since they do not lead to additional 
(algebraic independent) relations. 

Now, let us select some independent and recursive relations from the above 
ones, but first we take the basic identity (72) and the expression (54) for the 
speed in the spherical metric, whose zero order terms yield the two identities : 

Eo : = 46o, &o = ^ (97) 

E'o:vo = \j^l + 9oe\ (98) 

We shall use the symbols Ji , J2 etc. to denote various expressions which are of 
intrinsic type. By using (97) the identities Exq and £^40 can be restated as 

El : po(p? - 2/1) = Ji, Ji = 2uo - 60 (99) 

Next, the direction tpQ of 7* at the point [ifQ, 6q) is intrinsic; it is also conve- 
niently represented by the unit tangent vector 



* 1/ 9 ^ d . ^ d ^ d 

vq dip oO o(p oO 

The coefficients J^p, Jg are intrinsic functions, depending on the coordinate 
system, and they are related by the identity 



Therefore, we adjoin to our list (99) the two identities 



£2 : y'l = J^vo (101) 
E3 : 9i = Jevo 
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Still, wc have not used all zero order relations, namely £'20 and E30, and 
now we state them as 

E5 ■■ pI<P2 + PoPiVi = J5, Jb = 2/io + /o4®o (102) 

Eq : Po^2 + PoPi^i — "^6! J& = — 2—Ju,J0@o 

go 90 

By continuing this way, we obtain for each n > three new relations 

2 , 



E3n+2 : = p^<p„+2 + (103) 

Esn+S '■ = Po^n+2 + 

involving at each step the new triple p„_|_2, <y5„_|_2 , of variables of order n+2. 

Claim 4.2 It is possible to solve the above recursive relations for the variables 
(93) completely in terms of the intrinsic local geometric data in the shape space. 

This will be finally settled at the end of the subsection. At this point we 
have altogether 3n+8 variables 

P0> '^0; Pi, <^l, ^1; P2) V2) ^2; P„+2, <Pn+2> ^n+2; 

involved in 3n+8 recursive relations, and the first eight involve only the variables 
up to order 2. However, E[y E2 and are obviously algebraic dependent due 
to the identity (100), so let us search for one more independent relation among 
the variables of order < 2. Since we expect such a relation to involve local 
intrinsic quantities of order (at least) 3, a natural approach is to differentiate 
the basic identity (72) involving the Siegel function. Thus evaluation of the 
resulting identity (73) at f = to yields 

3^ + 2^ = Jr, Jr=^ = (^-%^) (104) 
PqVo 60 Wfe Afe 

Using the expression in (92) for vi we can restate the above identity as 



3^ 



V>if2 + -^fofiffi + goOi92 



J7 (105) 



PoVq 

By simple calculation and substitution using some of the previous relations Ei, 



pI 



1 2 

V>i^2 + ■^fo'Pi^i + 9o6\B2 



= 'fiAJb - PoPifi) + ^A\foJePlvl) + 6'i(s'o^6 - QopIpiOi) 

= -^pli^pl + goOl) + Vi{J5 + foJe&o) + OiQoJe 
Po 

= -4—60 + vo [J^Jn + foJ^Jg&o + goJeJe] 
Po 
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and by substitution into (105), using the identity PqUq 
sions for J5, Jg in (102), this leads to our new identity 



46o and the expres- 



Pi 

PqVo 



Js = 2eo\J^Ho + J0Vo)-J7= -^(2ui-6i) (106) 

too 



where ui is the tangential derivative of U* at Pq, cf. (90). 
Prom the system of algebraic equations 

Eq, E'q, Ei,E[, E2, E3, E4, 

we can now solve recursively and thus determine the variables 

Po^'"0,Pl,Vl,Ol,P2,V2,^2, 

successively in terms of the intrinsic data. In fact, this is obvious from the 
structure of the equations, once we have determined pQjVo, Pi, namely using the 
three equations Eo,Ei,E[: 



plvl = ieo, Po{pi-'2h) = Ji, 



Pi 



= Js 



It follows that 

2hp„ = - 

and consequently, for /i 7^ 0, 



- Ji + Jiplvl = - Ji + 4j|6o 



Po 



60 



(4J| + 1) - uo 



3/2 

Po 



&0, Pi 




In the case h = the identity (108) merely tells us that 

«o = ^(4J| + l)6o, 



(107) 
(108) 

(109) 
(110) 



and we can freely choose any initial size Pq of p and then calculate vq and pi from 
(109). In particular, we calculate the pair (<^j,^i) using (101), where the pair 
{Jtp, Je) represents the initial direction and hence is intrinsic. Now, we are able 
to calculate successively each new triple {p^,f^,6n), n = 2,3,4..., expressed in 
terms of the intrinsic data, as claimed above. 

Remark 4.3 The initial direction, {J^, Je), is the only basic intrinsic data with 

no invariant description, thai is, independent of the coordinate frame. However, 
from the recursive procedure it follows that Un,Vn, Pn,n > 0, come out with 
coordinate free expressions involving only the coefficients in (89). In fact, we 
can calculate in terms of p^ [i < n) and, Vj (7 < n — 1), by repeated 

differentiation of (72), next we calculate Un by applying differential operators 
such as (56) to U* , and finally pn+2 calculated using equation E\n- The 
beginning terms are 



Ui = VqUi, 



U2 



2- 1 

V0U2 + -VlUi, 



Vl 



46i - 3pgt;oPi 
2pl 



33 



4.2 Some basic results on the shape curves of three-body 
motions with vanishing angular momentum 

Let us first review some of tlie above facts from tlie local analysis and then 
draw a few immediate but important consequences. The above power series 
developments amount to the explicit calculation of the solution t 7(f) of 
the system (84) in the moduli space with the initial data (82). In doing so we 
started from the following 5-tuple 

(¥>o,^o),(^o,6o,ei), (111) 

where &q,&i may be replaced by Kq, Ki (or Kk, Kk+i, for the smallest k with 
Kk 0) which consists of three specific intrinsic local geometric invariants at 
the point {(pQ,9o) on the shape curve. In particular, we also recover the time 
parametrized shape curve 7*(i) by projecting 7(f) to the 2-sphcrc. 

Actually, since 7(t) is uniquely determined by the initial value problem, (82) 
and (84), it suffices to recover (82) from (111), namely the missing information 
in (82) is /Oq and pi. This turns out to be possible when /i 7^ 0, but for a "good" 
reason (see below) it is impossible when h = since in this case the shape 
curve only controls the product p^v^ = 46o. In any case, with the quantity 61 
we can actually determine p^ and, in particular, the question in Problem 3.20 
concerning the expansion index is settled. 

The general three-body problem has the 10 classical conservation laws (linear 
and angular momentum, and energy) due to its invariance under the Galilean 
symmetry group. All of them have been used and, in particular, the set of 
solutions is invariant under time translation, t t+to, as well as reversal of time 
{t —t) which reverses the direction of the trajectory. However, there is also 
an additional 1-parameter size/time scaling symmetry group, whose induced 
action on parametrized moduli curves sends 7(f) = {p{t),^*{t)) to 

7mW = iPir){t),l(r)it)) = {rp{r-h),r{r-n)), Vr > 0, (112) 

and changes the energy from h to = r~^h. In particular, although scaling 
and time translation leaves the oriented shape curve geometrically unchanged, 
its time parametrization is subject to an affine transformation 

t^at + to, a = r"^/^ > (113) 

Since the energy level /i = is scaling invariant, this also explains why the 
reconstruction of a unique initial size pg fails when h = 0. 

Remark 4.4 For any n > 1 the Newtonian n-body problem, has the above 1- 
parameter symmetry group {^r}, a,cting on size and time but leaves the shape 
invariant. For example, for a periodic motion with period P(i) and average (or 
initial) size the group sweeps out a periodic motion with the same shape, 
and the ratio P^^^/p^^-^ is independent of r. The case n = \ means the restricted 
case n = 2 with one of the masses (e.g. a planet) infinitesimal small, in which 
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case there is only one shape (a point) and the above ratio depends only on the 
large mass (the sun). This gives Kepler's third law, so the above symmetry group 
is essentially the generalization of this law. 

On the other hand, for three-body motions with vanishing angular momen- 
tum, the identity p^v'^ = 4(3 of (72) gives another quantity, p^v^, which is 
invariant under the above symmetry group. 

Definition 4.5 A time reparametrization 0/7(4) or 7*(t) by an affine trans- 
formation (113), for any a^Q, is called canonical, otherwise it is called excep- 
tional. 

Of course, in order to stay at a given nonzero energy level a canonical 
reparametrization must have a = ±1, and moreover, the orientation of the curve 
is reversed if a < 0. Now we can state the following basic unique parametrization 

theorem : 

Theorem 4.6 A three-body motion with zero angular momentum is, up to con- 
gruence and canonical reparametrization, uniquely determined by its oriented 
shape curve on the 2-sphere. In fact, it suffices to know the direction and the 
first two Siegel numbers ©0, Si o,t o,ny regular point on the shape curve. In 
particular, there are no exceptional reparametrizations. 

As shown before, the theorem still holds with (©o, 61) replaced by the cur- 
vature numbers {Kq,Ki) if the point is regular of order 0, and this is, indeed, 
the generic type of points. 

A curve on the 2-sphere which is the shape curve of a motion with total 
energy h can also be the shape curve for some motion with any other energy h' 
of the same sign as h. Indeed, we find the other motions by suitable canonical 
reparametrizations of the given motion, and by (109) it also follows that the 
sign of h (viewed as a number 0, ±1) is an intrinsic invariant at the shape space 
level. More precisely, we have the following quantitative measurement of the 
energy type : 

Theorem 4.7 Let 7* be a geometric curve on the 2-sphere, with the Siegel 
function & (with respect to U* , as usual), and consider the function 

2 

+ 1)-U* 

along the curve, where denotes the tangential derivative. Ifj* can be realized 
as the shape curve of a three-body motion with vanishing angular momentum, 
then the sign of A is constant along the curve (whenever A is defined), namely 
equal to the sign of the total energy h of the motion. 

Corollary 4.8 A given oriented (geometric) curve on the 2-sphere can be time 
parametrized in at most one way, up to canonical reparametrization, as the shape 
curve t 7*(0 of a three-body motion with zero angular momentum. Moreover, 
the sign of the total energy of such a motion is determined by the local relative 
geometry o/(7*,C/*) at a (regular) point. 



2 ^ 



Gdr^ 



6) 
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The above uniqueness property of time parametrization of geometric shape 
curves, the minimal amount of geometric information needed to determine the 
shape curve, and the monotonicity theorem which we shall discuss in Section 5, 
are our basic tools for the imderstanding of both the local and global picture of 
shape curves representing three-body motions with vanishing angular momen- 
tum. The monotonicity property tells us the m-latitude function is monotonic 
increasing or decreasing until the curve turns back somewhere in the opposite 
hemisphere. Thus the curve resembles an "oscillating motion" between the up- 
per and lower hemisphere which never stops, unless it ends at a triple collision 
or escapes to infinity. The curve crosses the equator circle transversely, or it 
goes to a binary collision and bounces back (via regularization) to the same 
hemisphere. In Section 7.5 we describe the problem of how to construct such a 
curve by linking together its maximal monotonic segments. 



5 The monotonicity theorem for shape curves 
5.1 A closer look at the gradient vector field of U* 

The analysis of trajectories, moduli curves or shape curves describing three-body 
motions depends, of course, ultimately on the function U* , whose behavior 
is largely reflected by the geometry of its gradient field. In this subsection 
some useful facts are established which are beyond those simpler statements 
concerning the critical or singular points of U* . 

We shall apply vector algebra in the Euclidean model for the moduli space, 
namely with M = M'"* as the Euclidean space (cf. Section 2) and vectors denoted 
by boldface letters. Thus the shape space M* — 5^(1) consists of unit vectors 
p = (x, y, z), IpI^ = x'^ + y'^ + z"^ , and p • q denotes the usual inner product. Set 

m = ^^mi, m = mim2mz, rhi = m2mz etc. and ^^TOi = 1 

3 

ki = 2^=^=, /?. =cos H—— )=sm i(^— ) 114 

Vl — rrii rrii + rrii rrii + rrii 

bi = (1,0,0), b2 = (cos^3,sin/33,0), ba = (cos/Jj, - sin/32, 0) 

where < /3j < tt is the angle between the binary collision points hj and 
bfe, {i, j, k} = {1, 2, 3}. The Newtonian shape potential function is the restric- 
tion of U to the above 2-sphere, 

U*=t^ = ±^ cf.(31) (115) 

where the mutual distances Sj = Vju are normalized to / = 1, and hence by a 
formula of Lagrange 

I = Y^misj = l (116) 



36 



The basic behavior of U* is, of course, given by the 8 special points on the 
2-sphere, namely 

bi,e3,b2,ei,b3,e2; Po,Po (117) 

where the first six are cyclically ordered (eastward) along the equator circle E* 
representing degenerate m-triangles. The are poles where U* tends to oo, 
and the Euler points ej are the saddle points. Finally, the remaining two are 
the minima, say po lies on the northern hemisphere and Pq is the symmetric 
(mirror) image with respect to the equator plane. 

We can also use the points bj to describe the gradient field, as follows. Let dx 
be an arbitrary infinitesimal vector perpendicular to p (i.e. rfx GTpS'^). Then, 
on the one hand 

VJ7*(p)-rfx = Z7*(p + rfx)-J7*(p) (mod|dx|^) (118) 

and on the other hand, 

U*{p + dx) - U*{p) = > ' fc, ( (2 - 2(p + dx) • b.)-* - (2 - 2p • b,)" 



rfx (mod |dx|^) (119) 




Set 



B = B(p)=5:^^ (120) 



i=l 

Then it follows from (118) and (119) that the gradient of the function U* is the 
orthogonal projection of the vector B to the tangent plane of the sphere at p, 
namely 

Vf/* = B - (B • p)p (121) 

The characterization of the critical points of U* is, of course, well known. 
However, with the following lemma we also like to establish the identity (122). 

Lemma 5.1 Let po and Pg represent the pair of equilateral m-triangles with 
1 = 1 and with opposite orientations. Then po and pg are the minima of U* , 
and moreover 

b, ■ (b, - Po) = b, . (b, - p^,) = Jl!lf!l]L^ (122) 

(1 — rriijm 

Proof. The determination of the critical points away from the equator follows 
readily by Lagrange's multiplier method in M with the constraint / = 1. As 
coordinates we can, for example, use the individual moments of inertia Ij, but 
the calculations arc simplest in terms of the mutual distances = Sk using 
(115) and (116). This shows the minimum (s?, 82,53) of U* (on any of the 
hemispheres) satisfies the following set of equations with a multiplier A 



1 

2R) 



0-3 -A, 1=1,2,3 (123) 
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and hence, by (116), all the sides are equal to l/V^ . Moreover, by (115) 

ki TTli /—— 



|Po-bi| 5° 



mjmkVm 



,2 Amjmk 



2-2po-b, = |po-b,r = — =^ (124) 

(1 — mi)m 

po • bi = 1 



(1 — mi)rh 
and this gives (122). ■ 

Remark 5.2 It is also straightforward to check the identities 

] rW = r^ 7772 = l^H^ - rrii) (125) 

|Po-bi| Ipo - b,| * 

^(1 - mi)hi = (126) 

i=l 

To simplify the notation below, let us write 

«,) = (I-iri, , = ^j^,cf.(122, (127) 

*i(i) = im^(l-TOi)Vi(0 (128) 

where ■(/'i is defined for t < Ci,i = 1,2,3, and observe that the derivative of 
is strictly positive. The following two lemmas will be useful. 

Lemma 5.3 The expression (120) of B can he restated as 

3 

B{p) = Y,MU)bi (129) 

where 

^, = b,-(p-Po) (130) 

Proof. By (125) 

■g _ kibi _ y-v ki / |po ~ b. 



i=ilp-bi| i=i|Po-bi| V|p-b. 



3 . 3 

7712 



= E^(i-"^4 b..(b.-p) J 



i=l 

rhi A / bj • (p - po) 

— — > (1 - mi) 1 - - — r 

4 V bi-(bi-po) 

^X^(l-m,)^,(^Jb,=^vi.,(^,)b, 

i=l i=l 



38 



Lemma 5.4 Let p be a unit vector different from po and pg. Then B • (p — pg) 

is strictly positive. 



Proof. By the mean value theorem, there exists < < 1 



and we recall that the derivative of is strictly positive. By (129) and (130) 



Note that B is, by definition, a linear combination of {bi, h^}, thus lying 
in the xy-plane and consequently 



where p (resp. po) is the orthogonal projection of p (rcsp. po) in the .xy-plane. 
Clearly, the geometric meaning of the positivity of the inner product in the 
lemma is that the angle between B and the vector from po to p is strictly less 
than n/2. 

5.2 The relative geometry between VC/* and the conjugate 
pair of co-eixial families of circles associated to {pojPo} 

In spherical geometry, circles are the simplest kind of curves and they are char- 
acterized by the constancy of their geodesic curvature. Associated to a given 
pair of points, such as the minima {po,Po} of U* , there are two co-axial fami- 
lies of circles, namely the family of circles passing through the two given points 
and its dual family consisting of those circles which are orthogonal to all circles 
of the former family. We shall denote the conjugate pair of coaxial families of 
circles associated to {po, Pq} by and J^'. 

To a given p G '5'^(1) other than Po,Po. let us denote the unique circle 
of !F (resp. !F') passing through p by Tp (resp. F^). In fact, Tp is simply the 
intersection of <S'^(1) and the plane spanned by the triple {po, Pq) p}- Therefore, 
the tangent line of Tp at p is the intersection of the tangent plane TpS^{l) and 
the above plane. In fact, the following pair of vectors 



at p such that T (resp. T') is tangent to Tp (resp. Tp). The vector T defines 
the southward direction, that is, away from po, whereas T' defines the eastward 
direction along V . 




i=l i=l 



B • (p - Po) = B • (p - Po) 
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For example, in the case of uniform mass distribution, {po, Po} = {N, S} are 
the north and south pole of the sphere, whose conjugate pair of coaxial families 
of circles are the usual longitude circles (or meridians) and latitude circles. In 
this case 

T = 2 [(p • k)p - k] , T' = -2p X k 

and these are positive multiples of the coordinate vectors ^ and ^ associated 
with spherical polar coordinates ((/?, 6) centered at the pole N. 

Proposition 5.5 The inner product between the "southward" vector T and the 
gradient vector \/U* at a point p on the sphere is 

T-VU* = T • B = [p-(po-p[,)] [B • (p - Po)] (133) 

In particular, on the northern (resp. southern) hemisphere the angle (3 between 
T and VU* is in the range < /? < 7r/2 (resp. 7r/2 < (3 <'jt). 

Proof. By (121), T-V?7* = T • B, and clearly B • Po = B • Pq, consequently 



T • B = B • {p X [p X (po - p'o)] + p X (po X p'o)} 

B p 1 B po p Po .^3^. 

B • (Po - Po) P • (Po - Po) B • Po p • p(, 

= (B • p)(p • (po - p[,)) + (B • Po)(p • (p; - po)) 

= [P • (Po - Po)] [B • (P - Po)] 

Finally, by Lemma 5.4 the second factor in (133) is always positive, whereas 
the first factor is positive on the northern hemisphere and vanishes precisely on 
the equator circle. ■ 



5.3 The monotone m-latitude theorem 

Spherical polar coordinates {(f, 9) centered at the north pole N parametrize, of 
course, the latitude and longitude (meridian) circles, which constitute the pair 
of coaxial families of circles associated to the pair {N, S*} of geometric centers 
of the two hemispheres. However, instead of using the colatitude let us rather 
parametrize the latitude circles by the latitude in radians, — 7r/2 < A < 7r/2, 
namely A = 7r/2 — </? and hence A is positive on the northern hemisphere. 

For equal masses the pair of minima {po, Pq} of U* happens to coincide with 
the pair {N, S}, but this does not hold for non-equal masses. However, there 
exists a unique Mobius transformation which maps po to N , Pq to S and the 
equator circle E* to itself. Such a Mobius transformation maps (resp. T') to 
the family of meridians (resp. latitude circles). 

Definition 5.6 For a given mass distribution {toi, m2, ma}, the m-modified 
latitude o/p € S^(l) is defined to be the latitude in radians of the image of p 
under the above Mobius transformation, and it is denoted by A(p). 
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For example, A(po) = 7r/2, A(po) = —it/2, and A(p) = if and only if p eE* . 
Moreover, A(p) = — A(p') for any pair {p,p'} representing similar m-triangles 
of opposite orientations. 

For a given (smooth) curve 7* {t) on the sphere we shall consider the associ- 
ated function 

A^.(t) = A(7*(t)) (135) 

which records the m-modified latitude along the curve. It turns out that for 

shape curves representing three-body motions with zero angular momentum this 
function has a remarkable monotonicity property. Namely, it oscillates between 
local maxima where it is positive and local minima where it is negative, and 
between two such extremals it is monotonic. The only exceptions arise when 
the function is a constant, as described by the following lemma. 

Lemma 5.7 // the m-latitude function is constant along 7*, then 7* is 
an exceptional shape curve which is either a single point or is confined to the 
equator circle (cf. Definition 3.12). 

Proof. Assume 7* is not a single point and is confined to an m-modificd latitude 
circle F' diff'erent from the equator. We choose a spherical polar coordinate 
system ((/?, 9) centered at the geometric center of F'; hence if is also constant 
along 7*. By equation (ii) of the system (34), [/* must be negative along F', that 
is, the gradient VC/* is pointing inward along the circle. However, F' encloses 
Po (or Pq) and Proposition 5.5 tells us that VJ7* is pointing outward, so this is 
a contradiction. ■ 

Now, let us assume 7* is not exceptional as in the above lemma. We shall 
state and prove the Monotone m-latitude theorem : 

Theorem 5.8 Let 7*(t), a < t < b, be a segment of the associated shape curve 
of a 3-body trajectory with vanishing angular momentum, and let A^* as in 
(135) be the function recording the m-modified latitude along 7* [a, 6]. Suppose 
that a < to < b is a critical point of A-y. (i.e. X'^.,{to) = 0) or is possibly a 
singularity. Then A-y*(to) niust be a local maximum (resp. minimum) when 
7* (to) lies on the northern (resp. southern) hemisphere. 

Proof. We may assume the point qo = 7* (to) on the m-modified latitude circle 
F' = T'^^ is strictly inside either the northern or southern hemisphere, since 7* 
crosses the equator transversely or it hits a binary collision point and bounces 
back into the same hemisphere (by regularization). Moreover, by using the 
reflectional symmetry which reverses orientation we may reduce the proof to 
the case that qo lies on the northern hemisphere. 

There are two cases to consider; either qo is a cusp, that is, the speed vq of 
7* vanishes, or qo is a regular point {vq ^ 0) and hence the curvature function 
of 7*(t) is smooth at to. 

If a cusp is encountered at qo, it cannot be the critical point po of U*. In 
general, the nonzero vector VJ7*(qo) actually gives the outgoing direction of 
the cusp, which by Proposition 5.5 is directed "southward" and hence A-y» is 
strictly increasing (resp. decreasing) when 7* approaches (resp. leaves) qo. 
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In the other case, 7* and F' are tangent to each other at qo. By reversal 
of time if necessary, we may assume that the velocity vector of 7* at qo points 
in the "eastward" direction of F', that is, the positive direction of F' as the 
oriented boundary of the circular cap containing po . Geometrically speaking, a 
local maximum of A-y. at to means exactly that the geodesic curvature Kq of 7* 
at qo is strictly less than that of F', which is a positive constant ko- Thus, it 
suffices to show Kq < fco, and we claim, in fact, that Kq < 0. 

Suppose to the contrary that Kq is at least positive, and recall Theorem 3.2 
and its identity (50), but scaled with M* as the unit sphere. Then, on the one 
hand 

(sin^ a)KQ > 0, 
and on the other hand, the positive normal vector to 7* is 

-T(qo) 



|T(qo)| 

and, by Proposition 5.5, evaluation at qo yields 

-^ln{U + h) = —-^VU*-iy* 
du* ^ ' U + h 



1 T 

-VU*-— > 



This implies that 



(U + h) |T| 



2(sin2 a)Ko - \n{U + /i) > 



and hence contradicts the identity (50). Consequently, Kq cannot be positive 
and, in particular, we conclude that A-y. (to) is a local maximum. ■ 

Corollary 5.9 Suppose that 7*(t) is the associated shape curve of a 3-body 
motion with vanishing angular momentum, without ever encountering a triple 
collision or escape to infinity. Then it contains infinitely many eclipse points 
(resp. local maxima and local minima for A-y ), and they occur at alternating 
sequences of times. 



6 The asymptotic behavior at a triple collision 

Three-body motions leading to a triple collision have vanishing angular mo- 
menum, and their moduli curves ^{t) are exactly those geodesic curves in 
{Mfi,dsf^) leading to the base point O = {p ~ 0) as the limit. Therefore, it 
is also natural to review and study their basic asymtotic properties at the triple 
collision. 

The classical works of Sundman and Siegel tell us that the triple collision 
is the only essential singularity of three-body motions, and their asymptotic 
theorem, briefly stated as Theorem 6.5 below, gives a qualitative description 
of the behavior at the singularity. We mention here some major works in the 
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classical literature which have contributed to the understanding of the collision 
motions, namely Sundnian[14] , [15], Levi-Civita[8] , Siegel[10], [11], Siegel-Moser 
[12], Wintner[16]. Unfortunately, the proofs one finds in the above literature are 
rather long and difficult, and thus it is worthwhile to provide simpler proofs, as 
well as improvements of their results, in the setting of kinematic geometry. 

The asymptotic theorem is, in fact, a direct consequence of the asymptotic 
estimates of / and its lowest derivatives / and /, and since p = \fl is the kine- 
matic distance to the triple collision (base) point, these lower order asymptotic 
estimates are needed somehow for any proof of the above theorem. On the other 
hand, the theorem can actually be regarded as the geometric interpretation of 
such estimates. 

Along the way we shall also give remarks on the works of Sundman and 

Siegcl, and in the final subsection we shall apply Wintner's idea of using a 
logarithmic time scale to deduce the asymptotic formulae for the time derivatives 
of / of any order, cf. Theorem 6.8. 



6.1 Ray solutions as a model for the asymptotic behavior 
at a triple collision 

We begin with some vector algebra in the Euclidean space M^^^© 1^(2)® ^(3) 
= of all triples 5 = (31,82,33), or rather in the subspace of m- triangles 
defined by ^ m^ai = 0, equipped with the Jacobi metric (19), cross product 
and exterior product 

5x5' = ^ micii X bi e M^ 5 A 5' = ^ ai A b^- e A^M^ 

where the standard basis vectors A eg S M^.^ A M^^.^ has length ^mirrij . Some 
useful relationships between these operations are expressed by 



I^A^T = det 



\5h5'\ 



5-5 5-5' 
S-S' S' -S' 

|2 



ISl"" IS'l"" - {S ■ S'f > \6x6' 



^ mf jai X bij^ + ^ mimjiiij (/z^^- > 0) 



In particular, for a motion 6{t) with velocity 6{t) and individual angular mo- 
menta flj, we deduce the relations 



6A6 



2IT- -P > \n\^ 
4 - ' ' 



(136) 



These can also be interpreted in terns of the splitting of kinetic energy 

. 2 



IP 



+ 



<5A<5 
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where T^, T^ot is due to radial motion, change of shape, or (rigid) rotational 
motion, respectively. In the case of planary motions, equality holds in (136) 
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since = jo 1^/21. In particular, equality holds if H = since in that case 
the motion is planary, e.g. by a simple geometric argument. 

Ray solutions provide, of course, the simplest examples of triple collision 
motions. Here each particle moves along a fixed line through the center of 
gravity (origin) and hence J7j = for each i, and also T" = T'^ = 0. By 
Remark 3.3, the shape of the ray is a critical point of U* on the 2-sphere M*, 
namely a Lagrange point or an Euler point. In other words, the motion is either 
a homothetic deformation of an equilateral triangle or a degenerate triangle 
which is an Euler configuration. Let ^, be the value of U* at the above critical 
point. Then the Lagrange-Jacobi equation (5) reads 

/ = 2/i^+4/i (137) 

and hence the only variable of the problem, I it), is the solution of a 1-dimensional 
Kepler problem. 

The equation (137) can be solved explicitly, but we seek the solutions with 
the (singular) initial condition 7(0) = 0. In the special case of /i = 0, 

m=Kt^, K=i^)i (138) 

and for general h there is a formula t = Fh{I) which can be inverted and, for 
example, this yields a series development of type 

oo 

i=l 

To facilitate our study of asymptotic estimates in general, let us introduce 
the commonly used notation 

/-g 4=> \im f/g = 1 asi^O (139) 

Then all ray solution have the same asymptotic behavior at i = 0, in the sense 
that their time derivatives of I{t) at i = yield the same asymptotic formulae, 
beginning with 

I{t)^Kt^, /(t) ~ J ~ ^ift-i, (140) 

and clearly the higher order asymtotic formulae follow the same pattern, namely 

^ = 0,1,2,3,... (141) 

For a general triple collision motion, it was first realized that f2 must vanish, 
hence also vanishes. On the other hand, although was found to be the 
dominating kinetic energy, T"^ cannot vanish for a non-radial motion and may 
perhaps tend to infinity at some lower order of magnitude. However, although 
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the general asymptotic behavior is certainly more involved due to the change of 
shape, it turns out that the estimates (141) still hold, by Theorem 6.8. Moreover, 
a general triple collision motion has one of the above simple ray solutions as its 
asymptotic limit, according to Theorem 6.5. 

Exact information on the limiting behavior of the shape is not really needed 
to derive the asymptotic formulae (141) for the motion in the radial direction. 
The theory of Sundman and Siegel establishes the formulae only up to fc = 2, 
namely the asymptotic estimates (140). To proceed from A; = 0, 1 to A: = 2 
Siegel introduced the following function g{t) and proved its crucial property 

g{t) = (8/T - /2)f-2/3 ^0 as t ^ (cf. Siegel[ll], Chap. Ill, §1) 

which in our setting can be reformulated as 
. 2 

6 AS I 

— — = 2pT'' = -p^v^ ^0 as t ^ (142) 

where v is the speed of the shape curve. From our viewpoint, we recognize the 
expression in (142) as the Siegel function & of the associated shape curve 5*{t) 
on the 2-sphere, cf. (72). 

6.2 The results of Sundman and Siegel 

The following basic fact on the vanishing of the angular momentum of three- 
body motions leading to triple collision had already been stated by Weierstrass 
when Sundman first proved the following classical statement at the beginning 
of the 20th century. 

Lemma 6.1 (Sundman) The angular momentum ri is necessarily zero for a 
triple collision motion. 

Proof. First, by translation and (possibly) reversal of time, we shall rather 
assume (in all Section 6) there is a triple explosion at t = 0. Using the Lagrange- 
Jacobi equation (5), it follows from 7 — > that 7 — > oo and hence 7 > for 
t S (0, to) and to suitably small. Sundman discovered and made use of the 
rightmost inequality in (136). Namely, for a given value |r2| > it is not so 
difficult to see that I{t) has a positive lower bound. 

However, wc shall proceed with a slightly different proof since (136) also 
involves the individual momenta rij, and this will enable us to prove a stronger 
version of the lemma (see Corollary 6.4). For that purpose we set 

C(t) = (5 A (5 ^ Co = inf C for t € (0,to) 



Now, multiplying the Lagrange-Jacobi equation by 7 gives the inequality 

77 = 7(2r + 2h) = ic + 2hi + ^i"^ > ic + 2hi 

I AI I 
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which by integration yields 



to to . 

j iidt > j jCdt + 2h{Io - I) 

t t 

\il > - P) > Co H^j-) + 2h{h - 1) 

2hl - 2hIo + hii 
ln(7o/7) 

Hence, Co = and the constant sum n = ^ f2 j is zero. ■ 

Remark 6.2 Clearly, the above proof also gives inf \ fti\ = for each i, but one 
cannot yet conclude that fli ^ 0. See Corollary 6.4 for this last step. 

The major results of Sundman and Siegel concerning a general triple collision 
motion d{t) can be summarized as follows. 

• Sundman : t~i6{t) tends to a limit whose shape is that of an equilateral 
triangle or an Euler configuration. 

• Siegel : The magnified or "big triangle" t~iS{t) approaches a fixed m- 
triangle Sq in the Euclidean configuration space M. 

Standard references for proofs of these statements are Siegel[ll] and Siegel- 

Moser[12]. One finds that results proved by Sundman are, typically, seen to 
express properties at the moduli space level, that is, statements about the mod- 
uli curve S{t). Siegel improved his results by lifting them up to the configuration 
space level, where he studied the motion of the "big triangle" in the Hamilto- 
nian setting and performed a series of successive canonical transformations to 
simplify the analysis. 

A major step was to establish the validity of the above asymptotic estimates 
(140), and with the following proposition we shall provide a proof of this - in 
the spirit of Sundman and Siegel. Moreover, for the sake of completeness, in 
the last subsection we shall extend the proof to the higher order asymptotic 
estimates, fe > 2 in (141), using ideas due to Wintner. 

Proposition 6.3 For a three-body motion 5{t) with I{t) as t ^ 0, the 
asymptotic estimates in (I40) hold, where K is the expression in (138) with 
fj, = \imU*{5*{t)). 

Proof. For a ray solution, I = ^ Kti and pp^ ~ ^Ki . This suggests a 

study of the asymptotic behavior of pp^ for triple collisions in general, using 
equation (5) and kinematic geometry. Now, = and the kinetic energy has 
the splitting 

T = lif/ = T^ + T^ = + 1^2^2^ (25) (143) 
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and we set 

m{t) = pTp, 6(0 = ^e(t) = pr- (i44) 

where & is the Siegel function, see (63), (72). By the energy integral U* = 
pT — ph we can also write 

U* ^ 1R{t) + 6{t) as t ^ (145) 

Our first claim is that 

^R{t) p>0 and 6(t) ^ (146) 

as t ^ 0. Let us differentiate $H and substitute for I using (5), or equivalently 
(i) in (34), to obtain 

[°mt = 9l(io) - fH(i) = / ° pT^dt + h{p{to) - p{t)) (147) 
Jt Jt 

Since 91 > and the integral on the right side is > 0, linifH(t) = /u > must 

exist, that is, d\ = p + o(l) for small t. 

Suppose we had p — 0. Since min U* > 0, equation (i) in (34) implies 

pp2 = [/* + 2hp -29\=U*+ o(l) > C > 

and consequently 

PP=-^{pp') = ^{U*+o{l))>C-^ 

p{hf - p{tf =2 jj ppdt >'^C jj ^dt 

p{t) p{to)' 

This is clearly impossible, so we conclude p> 0. 
Next, we deduce successively 

ipp2^^ + o(l) ^ p^ = h2p + o{l))^ 
z p 

p = ^ p-^ +o{p-^) (148) 



/ y/pdp= I 

Jo Jo 



Py/pdt 



{^J2p + o{l))dt = ^2pt + o{t) 



(149) 
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and hence by (148) and (149) 



(150) 



6 Z p 



(151) 



Next we show &{t) 0. The integral on the right side of (147) exists, but 
the integrand is 

^ p 3 t 

and the integral of 1/t is divergent, hence liminf (S(t) = 0. 

It is, however, more difficult to show limsup(3(t) = 0, but let us apply an 
idea from Siegel[ll]. Namely, suppose to the contrary, that limsupS(t) > 0. 
Then, for a given e > there is an infinite decreasing sequence of numbers in 
(0, to), e > ti > t2 > ■■ > tk > 0, limtj = 0, so that 



e < &{t) < 3e, for t e Jk = [t2k, t2k-i] 
S(i2fe)=e, 6(i2fc-i)=3e 
\Hhk)-Ht2k-i)\ <e 



(152) 



By (145), in each interval Jk, U* and hence also the norm of Vf7* are bounded 
by the same constant C\ , and then it is not diflacult to show 

t = u = o{ri), forie Jk 

But for small t we also have 



T = U + h = -{U* + hp)< — 
P P 



for a suitable constant C2, consequently 

pf + pT 



C 

< —, for t€Jk 



(153) 



(154) 



By (152) and (154), 

2e = 6(t2fe-i) - &{t2k) = {pT - ^)t=t,,., - {pT - m)t 



< Co 



f*^"-^ dt 
Jt2k * 



+ e 



which implies 



rt2 



t2k-l /■*2fc-l 

^dt > e / 



— > ^ , for each k 
t C2 



(155) 
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and hence the sum of the integrals is infinite. On the other hand, 



6{t) pT" 



-pT" 



and the integral of pT"' on [0, to] exists by (147), so this is a contradiction. 

Having proved that &{t) 0, it follows from (145) that U* p, and now 
the Lagrange-Jacobi equation yields 



4 2 



I r^2T '^2TP = p^ 
where K is the expression from (150). This completes the proof. 



(156) 



Corollary 6.4 The quantity 6 Ad in (136) tends to zero at the triple colli- 
sion. In particular, the individual angular momenta fij as well as the "mixed" 
momentum term ^mix tend to zero. 

Since &{t) 0, the above statement follows immediately from 



SaS 



2IT - -P = 2IT'' = 2pe 0, cf. (136) (157) 



By "infinite magnification" at the triple collision the solution 5{t) coincides 
with one of the ray solutions in Section 6.1. This is the idea behind the classical 
asymptotic theorem, and now we give a simple proof of this in the setting of 
kinematic geometry. 

Theorem 6.5 (Sundman-Siegel) Any triple collision orbit is asymptotic to one 
of the ray solutions. 



Proof. By (51) and (156) 



2T: 



^dt' 



cos^ a 



4 

^-Kt- 

(dp 2 
^ds' 



(dp 2 
^dt' 



2TP 



1 as t ^ 



(158) 



which simply means that the moduli curve 6{t) of the given triple collision mo- 
tion 5{t) is tangent to a ray or, equivalently, the limit of its infinite magnification 
exists. It also follows that the limit ray must itself be a geodesic in Af/i, namely 
one of those rays representing the shape of a Lagrange or Euler configuration. 

However, the claim is also that 5{t) itself approaches a ray in M. To see 
this, consider as above the angle a between the radial and tangential direction 
in M, that is, the angle between the vectors 5{t) and 5{t). It follows that 



5-5 



cos a 



1^1 



2T 



1 as t ^ 



(159) 
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The limit ray in M projects to a geodesic ray in Mh, namely a ray consisting 
of the homothetic images of either a fixed equilateral triangle or a fixed degen- 
erate triangle of Euler's type. As shown in Section 6.1, these are the rays which 
admit triple collision motions, and thus the given motion S{t) will be asymptotic 
to the corresponding limit ray solution with the same energy h. ■ 

Corollary 6.6 The shape curve S*{t) converges to a Lagrange or Euler point 
5q on the 2-sphere, and the "big triangle" t~'^/^6{t) converges to an m-triangle 
Sq with the shape Sq and moment of inertia 

io = K = {^)i,^ = U*{S*) (160) 

Remark 6.7 Actually, a limiting shape of Euler's type cannot be reached unless 
the whole three-body motion itself is collinear, see e.g. §i5 in Siegcl-Moser[12], 
which refines and improves the classical Sundman-Siegel approach. The latter 
is described in detail in Siegel's lectures [11] of about 240 pages. On p. 138 he 
writes :"The difficulty of the problem consists in the fact that we cannot yet 
prove (this will be proved only at the end) that the big triangle referred to a fixed 
coordinate system has a limiting position as t ^ 0; all that we have proved so 
far is the existence of a limiting configuration relative to a rotating coordinate 
system. The triangle itself may go on rotating about its centre of gravity, ..." 

However, although finiteness of the rotation of the "big triangle" was proved 
via the convergence of the "big triangle" , neither an estimate of the actual angle 
of rotation nor its precise definition was addressed in the above studies of the 
triple collision. In reality, the "big triangle" is approaching its final shape and 
position quite fast and in a monotonic way. To make this precise, we propose to 
measure how much the equilateral limiting triangle So deviates in position from 
some natural reference equilateral m-triangle C, depending on the given collision 
motion 6{t), but also C = C{t) will be a function depending on the chosen time 
interval [0, t] under consideration. This goes as follows. 

We may assume the shape curve S*(t) is on the northern hemisphere and 
hence starts at the Lagrange point po = Sq. Let Si = S{ti) be the m-triangle 
at a given time ti > and write SI = S*{ti). Then there is a unique linear 
m-triangle motion 

z{t) = ^^^Ci + ^Si, te[0,ti] 
tl tl 

with vanishing angular momentum, connecting Si to some equilateral m-triangle 
Ci = C(^i) (cf- [5] , Section 3.3). The shape curve of this (virtual) motion is 
the geodesic arc on the sphere 5^(1) from the Lagrange point po to the point 
Si , and together with the curve segment of S* (t) from S^ to po they constitute 
a closed curve Ci on the sphere. We define the rotation angle tp{ti) of S{t) at 
the triple collision, measured from time t = ti, to he half of the signed area 

V'(ii) = ^Area{Di) = [ to (161) 
^ JCi 
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of the region Di enclosed by Ci. This is motivated by the kinematic Gauss- 
Bonnet theorem (of. [5]) for three-body motions with zero angular momentum, 
where traversal of a loop on the 2-sphere amounts to a net rotation (i.e. a 
geometric phase ) of the m-triangle in the configuration space, which can be 
calculated as the line integral of a kinematic 1-form uj (depending on the mass 
distribution and region of 5^). Moreover, duj = \dA where dA is the area form 
of the unit sphere. 

Now, for ti not too large, the shape curve 6*{t) will stay on one side of the 
geodesic arc since its curvature will have a fixed sign. So the rotation angle (161) 
decreases monotonically to zero as ti 0, and hence the kinematic geometric 
approach explains Siegel's angle of rotation and yields as well a recipe for how 
to measure it quantitatively. 



6.3 Higher order asymptotic estimates at a triple collision 

The asymptotic formulae for the energy functions S = T, T^, U and their time 
derivatives up to order k can be developed inductively together with those for- 
mulae for I up to order k + 2. However, from the three identities 

7"=^, T = U + h = \i-h (162) 

it is easy to show that the three cases of S, for a given order k, yield the same 
asymptotic formula. Therefore, the final description of the asymptotic behavior 
of the above quantities can be stated as follows. 

Theorem 6.8 For a three-body motion S{t),t > 0, with a triple collision at 
t = 0, the following asymptotic estimates hold as t ^ : 

where K = , fJ, = U*{6q), and 6^ is the limiting shape. 

From the initial cases fc = 0, 1, 2, proved in the previous subsection, wc shall 
complete the proof of the above theorem for A: > 2 by deducing the following 
equivalent formulae for the behavior of p = VI, 

^p^^(V^P), .>3, (163) 

In fact, they will follow inductively as a rather direct consequence of Newton's 
equation (2) and its energy integral, namely 

VU{6) = -^S, T = U + h, (164) 

but only after an appropriate transformation of space and time. This is the 
composition of a time dependent space transformation and a pure time trans- 
formation, as follows : 
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• Magnification of the motion 5{t) = (ai(t), a2(t), a3(t)) by the time factor 
t"^/^, as in the works of Sundman and Siegel, to assure convergence at 
t = of the magnified motion. We use the notation 

~5{t) = t-2/3<5(i) = (ai 52, ag), a,(t) = f-^^^a^it) (165) 

m = mt)), m = fm), 5o = iim<5w aee) 

where / is any (homogeneous) function on M or its tangent bundle which 
we shall evaluate along the trajectory. 

• A logarithmic transformation of time; set 

u = —logt (oT t = e^") 

and hence t ^ means u ^ oo. This transforms a function g{t) to the 
function g{u) = g{e~'"). 

The composition of the two transformations yields the motion u — > 6{u) in 
M, and we write 

d{u) = ^(e-") = (ai,a2,a3), aaiu) = ei«ai(e-") (167) 

f{u) = fiHu)), p{u) = \s{u)\ = ei«p(e-«), = 

This motion is, of course, a solution of the transformed equations in (164), which 
can be stated as 

with the appropriate interpretation of T, U and VU, cf. (172), (173). For 
example, from the above definitions 

U{u) = U{S{u)) = e-i«J7((5(e-")) = t^U{t) 

To derive the above equations and prepare for its usage, we shall make a few 

more definitions and establish some useful identities for differential operators 
generated by ^ and For functions of t (or u) it is convenient to write 

/i ~ /2 ifi - 12) = oil) as i -> (or M ^ oo) 
and, for example, since the magnified motion converges, 

5k5q, p« |5o| = ^^K (170) 

We say f{t) has order g at f = if 
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and then the notation 



m = t-^f{t) « /o 



(171) 



is consistent with (166) since a homogeneous function g of degree d on M, with 
5(^0) 7^ 0, has order q = 2rf/3 at f = 0. 

The transformed potential function, kinetic energy, and gradient are given 

by 



IJ = U{5) = Y^ 



niimj 



i<j 



du 



\ nil oa.1 I 



(172) 
(173) 



and by substituting these expressions together with 



2-^ I d 



du^ 



into the equations (164) one obtains the system (168)-(169). 

Lemma 6.9 // f{t) has order q at t = 0, with f{t) ~ /o and q ^ {0, 1,2, ..}, 

then there is the equivalence 



l<k<m, t''^f{t)-Q 

\<k<m, §,m-§,iM^) 



Proof. By applying the Leibniz formula 



d"" 
dt" 



\ nn-i Ji 



one proves the above equivalence by induction on m. We refer to Lemma 6.1 in 
[13] for a detailed proof. ■ 

Moreover, using the operator identity 



d^ d d'' 



associated with the logarithmic time change, t 
following equivalence 



(174) 

log t, one can verify the 



l<k<m, t'^—rg « 
- - dt''^ 



d^ 



1 < k < m, ——rQ ~ 



(175) 



The reason for introducing the change of variable f — > u is the following 
useful lemma of Tauberian type. 



53 



Lemma 6.10 (cf. #363 in Wintner [16]) Let f{u) he defined for u > Q and 
assume f{u) has a limit and ^/(w) is hounded as oo. Then ^/(u) 
as CO. 

Finally, we turn to the proof of the asymptotic formulae (163). The "initial" 
data needed to start up are provided by Proposition 6.3 and (170), which by 
Lemma 6.9 and (175) can be restated as 

p^^.ij^J^,±i,^^^p^, (176) 

Then, by equation (169), we first deduce T w 0, or equivalently w for 
each i. Moreover, by (176) each is bounded, tj is bounded and hence |ai — aj| 
has a lower bound when i ^ j. It follows that all partial derivatives of C/, with 
respect to components of an and of any order, are bounded (as functions of u). 
In particular, in equation (168) WU is bounded and hence also -^5 is bounded. 
Now, apply the operator ^ repeatedly to the equation (168) and deduce 

d'' - d^^"^ - 

——rVU is bounded => , , , .. S is bounded, for all > 1 

By the above Tauberian lemma it follows that 

d^ 

—rl w 0, for all > 1 

du'^ 

Similarly, apply ^ successively to the equation (169) and deduce that the high- 
est derivative of p is always bounded, hence by the Tauberian lemma 

—rp w 0, for all fc > 3 (177) 

aw* 

By the equivalence (175) and Lemma 6.9, the statement (177) is equivalent to 
the statement of (163), and this completes the proof of Theorem 6.8. 

Remark 6.11 The asymptotic estimates in Theorem 6.8 are also valid for a 

general collision (i.e. total collapse) of an n-hody motion, for any n > 2. The 
proof is esssentially the same as ahove and the previous suhsection, since we 
have used only the Riemannian cone structure of the moduli space M and, for 
exam,ple, the angle a is similarly defined for any n > 2. In fact, the actual 
structure of the shape space is irrelevant as far as the asymptotic hehavior of 
the radial motion is concerned. Moreover, the exponent u of t in the formula 
I ^ Kf^ is independent of n, hut depends on the degree — e of homogenity of the 
potential function, C/ ~ namely v = 4/(2 + e) where we assume < e < 2, 
and e = 1 is the Newtonian case. We refer to [13]. 
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7 A brief discussion of some open problems 



In this concluding section we shall formulate and explain some natural open 
problems in the present geometric setting. Recall that the trajectories of 3-body 
motions with zero angular momentum are already uniquely determined up to 
congruence by their associated moduli curves, which can be characterized (ge- 
ometrically) as geodesic cmwes in {Mii,dsf^). Furthermore, these geodesies to- 
gether with their time evolution are essentially determined by their shape curve 
on the 2-sphere. Finally, we recall two major results concerning shape curves, 
namely the unique parametrization theorem and the monotone m-latitude the- 
orem (cf. Section 4.2 and 5.3). Therefore, from now on shape curve means 
geometric shape curve, unless otherwise specified, and we also assume they are 
oriented. 

The geometric behavior of these spherical curves raises many interesting 
questions for an in depth understanding of 3-body motions. Here we propose a 
few natural problems of basic importance. 

7.1 Shape curves of periodic motions with vanishing an- 
gular momentum 

The study of periodic orbits is naturally a central topic of the 3-body problem 
as a whole. Clearly, the moduli curve and the shape curve of such a motion 
are periodic, but the converse may not be true. Therefore, we say the three- 
body motion is congruence periodic or shape periodic if the moduli curve or 
shape curve, respectively, is periodic with respect to time. However, it is an 
important consequence of the unique parametrization property that the time 
parametrization of these curves is dictated by the geometry of the shape curve, 
whenever the latter is non-exceptional. Then the notion of congruence periodic 
is the same as shape periodic, and this means the shape curve is periodic in a 
geometric sense which we explain as follows. 

Since it is natural to allow binary collision points, periodic shape curves 
can be characterized as the topologically closed shape curves. Namely, the 
curve is either the immersion of a circle, and we call it circular periodic, or the 
immersion of a closed interval (of length > 0) and is contractible, and we call 
it string periodic. In the latter case the curve is a "string" with two end points 
which are either a reversing cusp (i.e. at the Hills's boundary) or a collision 
point, and in order to qualify as a periodic curve it is tacitly assumed we take 
two copies of the "string" with the opposite orientation. 

Remark 7.1 A shape curve consisting of a single point p is an exceptional 
case, and there are additional string periodic moduli curves of the fixed shape 
p. Namely, p must be a Lagrange or an Euler point, and the ray solution with 
negative energy, starting at rest from {po,p) on the Hill's boundary (cf. Section 
7.4), leads directly to the triple collision point O. By traversing this ray segment 
in both directions we obtain a string periodic moduli curve. 
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Next, we shall describe the distinction between periodic three-body motions 
and shape periodic motions (i.e. of circular or string type). Let 7* be a closed 
(piecewise smooth) curve on the 2-sphere which is the shape curve of a motion 
7(^)1 to < t < ti, of m-triangles with vanishing angular momentum, and assume 
7* is periodic as above. Then j{to) and 7(fi) are congruent m-triangles and 
hence differ only by a rotation angle A^, which we can calculate as a line integral 
along the shape curve, according to the kinematic Gauss-Bonnet theorem, see 
(161). 

In particular, Aip is zero if the shape curve is string periodic, and hence 
the given motion "/{t) must also be periodic. Thus the notions of "periodic" 
and "shape periodic" are identical in this case. On the other hand, a circular 
periodic curve 7* encloses a signed area AA (depending on orientation and self- 
intersections) , and the above line integral can also be expressed as a surface 
integral which yields 

At/j = ^AA (178) 

Therefore, the motion j{t) is periodic if and only if the angle (178) is a rational 
multiple of 27r, say Aip = {p/q)2'jT with {p,q) = 1, and hence the number 
<]{ti — to) is the period of the motion. 

Thus the study of periodic 3-body trajectories is completely reduced to the 
study of closed shape curves on the 2-sphere, and it is a challenge to describe 
or characterize the various types of these curves in terms of simple geometric 
invariants. For example, due to the monotonicity theorem it is natural to regard 
the number of eclipse points (counted with multiplicity) as a measure of the 
complexity of the curve, hence the simplest curves are characterized by a small 
number of eclipse points. 

Problem 7.2 What is the minimal number of eclipse points on a (string or 
circular) periodic shape curve ? What even numbers can be realized? What are 
those periodic curves with a small number of eclipse points, say up to 10? 

The homotopy classes of closed cm-ves inside P = — {bi, b^} arc ele- 
ments of the fmidamental group 7ri(P), namely the free group of two generators. 

Problem 7.3 What are those homotopy classes of closed curves in P which 
can be represented by circular periodic shape curves ? 

Definition 7.4 We propose to define the chaoticity of 7* to be the following 
value 

' Area{S^) 

where -D(7*) C S'^ is the closure of the set 7*. and we say 7* is chaotic or 
non-chaotic if ch('-f*) > or ch('-f*) = 0, respectively. 

Problem 7.5 What are the possible values of chaoticity for shape curves rep- 
resenting motions with CI = ? 

Problem 7.6 What are the non-chaotic shape curves other than the periodic 
ones ? 
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7.2 Triple collisions 



The works of Sundman and Sicgcl show that triple coUisions is the only type 
of essential singularity of 3-body motions, while the binary collisions can be 
regularized analytically [15], [8]. Recall that the singularities of the Newtonian 
potential function U in the moduli space M ~ consists of the triple of 

rays |Obj, i = 1, 2, 3|, where the base point (or origin) O represents the triple 

collision and the rays represent the three types of binary collisions. Moreover, 
the moduli curves of triple collision motions with total energy h are exactly 

those geodesic curves in {Mh,dsfJ with the point O as a limit. The following 
arc some pertinent problems on the geometry of such geodesies. 

Problem 7.7 The existence (resp. uniqueness) problem on the shortest path in 
{Mh, ds\) linking a given point p in Mh to the base point O. 

By the scaling symmetry we may assume h = 0, ±1, where the case of 
/i = — 1 is most difficult and also most interesting. The existence of such a 
shortest geodesic curve between O and a given point p in Mh can be proved by 
Hubert's direct method when /i = or 1, whereas for the case oi h = —1 the 
existence will depend on the position of p in the Hill's region Mh- Note that 
Hubert's direct method also applies for those p with 

dip, O) < dip, dMh) + didMh, O) 

The uniqueness problem is, however, much more interesting than the exis- 
tence problem, but it is also much more difficult and subtle. For a geodesic 7 
starting from O, the question is how far out 7 is the unique shortest geodesic 
from O. We remark that for points lying in the eclipse plane, the shortest 
geodesic is not in the eclipse plane, and hence the limiting shape at O of the 
shortest geodesic must be a Lagrange point, say po. By the monotonicity of its 
shape curve 7 will eventually reach the eclipse plane, but after the first eclipse 
7 ceases to be of shortest length. Hence, the best we can hope for is uniqueness 
up to the first eclipse point. 

We propose to investigate first the case of /i = 0, due to the scaling invariance 
of this energy level. Then the general uniqueness of a shortest geodesic between 
any point p and O reduces to the uniqueness for eclipse points p lying at the 
distance p = 1 from O, namely for points on the eclipse circle E* . Thus the 
problem is reduced from the moduli space M to the shape space M*, namely 
we ask about the uniqueness of such shape curves between points on E* and 
Po- 

Problem 7.8 For the case of energy level h = and for a given mass distribu- 
tion, let S be the set of triple collision moduli curves emanating from O, whose 
shape curve starts out from the Lagrange point po on the upper hemisphere of 
M* = S"^. Let S* be the initial arcs of the shape curves from po (but po not 
included) to their first point on the equator circle E* . Is the family of curves 
S* a foliation of the punctured upper hemisphere 5^ — {po} ? 
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In the case that there exists a unique shortest geodesic in {Mh, dsfj Unking 
a given m-triangle S to the collapsed m-triangle O, it is certainly interesting 
to actually estimate its length, initial direction and total rotation angle of the 
triangle in terms of its geometric invariants. See also the last part of Section 
6.2. 



7.3 Binary collisions and nearby trajectories 

The base point O is, of course, the only singularity for (M,ds'^). But the Rie- 
mannian manifold (Mh , dsf^ ) has another kind of singularity along the triple of 

rays |obi,7' = 1,2.3| mirms the initial point O, say, of binary collision type. 
One expects that understanding of the geometry of geodesic curves in the vicin- 
ity of this type of singularity will be an important topic in the study of the 
global geometry of geodesies on (M^, ds^). 

Problem 7.9 What kind of local analysis will enable us to provide an effective 
control on the local geometry of geodesic curves in the vicinity of a singular ray 
of a given binary collision type ? 

Let us make some further remarks. In the vicinity of the ray Obi the gradient 
vector field VU is closely approximated by the field Wi, where 

rjk \/l - mi di 

and di = di{x) = p sin cTi is the distance in {M,ds^) between x and the ray 
Obj, cf. (31). Thus, there is a suitable rotationally symmetric metric which 
provides a good approximation of ds\ when we are close to such a ray singularity. 
Application of Noether's theorem to this simpler metric yields a first integral 
of its geodesic equation which is almost constant along a geodesic segment of 
{Mh, ds\) near Ohi. This will serve as a useful auxiliary function whose analysis 
will provide an effective control on the above local geometry. 



7.4 Trajectories stEirting at the boundary of the Hill's re- 
gion 

In the case of negative energy, say h = —1, the variety Adh is the Hill's region, 
namely the proper subset of the moduli space from which the moduli curves of 
the three-body motions cannot leave. The region is enclosed by its boundary, 
namely Hill's surface dMh which is the smooth surface defined by p = U*{(p, 6), 
with the (kinematic) gradient field VU as a normal field. 

Here we shall focus on those geodesies of the metric dsf^ starting at the 
surface dMh, where the metric becomes identically zero. Hence, a curve lying 
on dMh has zero length, and a minimizing curve containing a segment on dMh 
is only virtual and cannot, of course, represent an actual trajectory of a 3-body 
motion. Therefore, in the study of variational problems of this kind one often 
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needs a certain estimate or geometrical control of those geodesic curves starting 
at dMh, that is, the moduli curves of those 3-body motions with no kinetic 
energy at t = to- They constitute a family {73;^} of geodesies parametrized by 
their initial points xq G dMh- 

Following Jacobi, it is natural to study the variational vector fields along 7^^^ 
with respect to variations within the above family of geodesies. These vector 
fields are solutions of the Jacobi equation along 7^^^ with their initial vectors 
belonging to T^„(dMh). 

Problem 7.10 Let .tq be a generic point of dMh, h — ~1, and let 7^^ be the 
geodesic curve with initial point xq. How do we obtain an effective (that is, 
simple and useful) lower bound estimate of the distance between Xq and the first 
zero point of Jacobi vector fields of the above type, in terms of the geometric 
invariants at xq ? 

7.5 On the problem of fundamental segments 

For a fixed energy level h = 0, ±1, consider the family S(/i) of all oriented geo- 
metric shape curves, with the exceptional ones removed, of three-body motions 
with zero angular momentum. According to the monotone m-latitude theorem 
the curve 7* can be viewed as a union of its segments Ci = {pi,qi), between 
two consecutive points Pi and Qi of extremal m-latitude. We shall refer to them 
as the fundamental segments. Thus the end points Pi,qi lie on opposite hemi- 
spheres, unless one of them is a binary collision point (and hence lies on the 
equator circle), and moreover, the m-latitude is strictly monotonic along the 
segment. Clearly, a global shape curve can be regarded as being pieced together 
by such fundamental segments, and a periodic shape curve has only a finite 
number of them. 

Conversely, we may try to construct curves by connecting d to Ci+i in a 
"smooth" way. Here Ci and C^+i belong to so "smooth" means their 

union also belongs to S(/i). For simplicity, assume we are using only regular 
fundamental segments C = {p, q) , that is, p and q are regular points. Ob- 
serve that C is tangential to the m-latitude circle at p, so its direction will be 
completely specified by an index £ = 0, 1 representing "eastward" or "westward" 
respectively. Thus we can associate to the starting point p = {'PtO) the following 
5-tuple of numbers 

[p] = (^,(?, 60,61, £) (179) 

which determines C completely and therefore also the 5-tuple [q] associated to 
its end point. In (179) 60, 61 are the Siegel numbers of C at p, as explained in 
Section 4.1 and 4.2. 

Roughly speaking, the relationship between the initial data and terminal 
data for a fundamental segment with regular end points provides a type of 
correspondence 

[p] = (<p, e, 60, 61, £) ^ [q] = (tt - e', 6'o, e[,s') 
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on a dense open set of S"^ x M+ x Rx {0, 1}. Moreover, since p and q arc points 
on opposite hemispheres, let us compose the above correspondence with the 
reflectional symmetry with respect to the equator circle, namely we replace [q] 
by [p'] = (lp\ 9' , 6q, &[,e'). Finally, wc assume p (and hence also p') lies on the 
upper hemisphere, thus arriving at the fundamental correspondence 

eiToUTi^ToUTi, (180) 

where the 1i are identical copies of the 4-dimcnsional space T = S"^ x ]R+ x M. 
The correspondence is defined on a dense, open set, where it is also invertible. 
In fact, with some more labour it would be possible to extend the fundamental 
correspondence to include irregular points (i.e. cusps and collisions) as well. 

Remark 7.11 A periodic shape curve is the assemblage of a finite number of 

fundamental segments whose initial data constitute a periodic orbit of the above 
correspondence (180). Namely, if [p] has even order 2k, then the orbit of [p] 
defines 2k fundamental segments which join together to a periodic curve. On 
the other hand, if the order is 2fc + 1, then the end of the curve lies in the 
southern hemisphere, so by running through the orbit twice the order will be 
Ak + 2, and the associated curve will be periodic. 

Thus, the correspondence (180) provides a natural way to a systematic study 
of the geometry of global shape curves. 
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